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ABSTRACT 

In this paper the current release of the Monte Carlo event generator Sherpa, version 1.1, is presented. 
Sherpa is a general-purpose tool for the simulation of particle collisions at high-energy colliders. It 
contains a very flexible tree- level matrix- element generator for the calculation of hard scattering processes 
within the Standard Model and various new physics models. The emission of additional QCD partons 
off the initial and final states is described through a parton-shower model. To consistently combine 
multi-parton matrix elements with the QCD parton cascades the approach of Catani, Krauss, Kuhn and 
Webber is employed. A simple model of multiple interactions is used to account for underlying events 
in hadron-hadron collisions. The fragmentation of partons into primary hadrons is described using a 
phenomenological cluster-hadronisation model. A comprehensive library for simulating tau-lepton and 
hadron decays is provided. Where available form-factor models and matrix elements are used, allowing 
for the inclusion of spin correlations; effects of virtual and real QED corrections are included using the 
approach of Yennie, Frautschi and Suura. 
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Introduction 

The LHC will pose new challenges to both the experimental and the theoretical community. It will 
operate at the highest centre-of-mass energies ever reached in a collider experiment and provide an 
enormous luminosity, leading to tremendously large event rates. On the experimental side, the huge 
phase space in conjunction with the exciting and intricate physics programme of the LHC, ranging from 
high-precision flavour physics at comparably low scales to the discovery of new particles in the TeV 
range necessitated the development and refinement of triggers and analysis techniques. In addition, 
from a more technological point of view, data acquisition, storage and processing therefore required the 
creation of a world-wide network satisfying the greatly increased computing needs. On the theoretical 
side, on the other hand, demands for higher precision to correctly model signals of new physics and their 
backgrounds led to the rethinking of calculation and simulation paradigms and to the construction of 
a new generation of modern tools. Maybe the most prominent manifestations of these paradigm shifts 
are the new, full-fiedged event generators, which certainly will prove to be indispensable tools for data 
analysis. Currently, these new simulation programs holding many new features are replacing the well- 
established traditional ones. In many cases, they allow a wider range of applications; for example, the 
typically more modular frameworks have alleviated the incorporation of new physics models. Often the 
new tools also offer higher precision in the simulation, because better and more accurate techniques at 
various stages of the event generation have become available over the past years. Ultimately, this led to 
a drastic improvement e.g. in the description of Standard Model backgrounds to signals for new physics. 
Thus it is no surprise that there is a widespread belief that these new tools will have a sizable impact 
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Fig. 1 Pictorial representation of a tih event as produced by an event generator. The hard interaction (big 
red blob) is followed by the decay of both top quarks and the Higgs boson (small red blobs). Additional 
hard QCD radiation is produced (red) and a secondary interaction takes place (purple blob) before 
the final-state partons hadronise (light green blobs) and hadrons decay (dark green blobs). Photon 
radiation occurs at any stage (yellow). 



on the understanding of LHC physics. The construction, maintenance, validation and extension of event 
generators is therefore one of the principal tasks of particle-physics phenomenology today. 

The inner working of event generators 

Fig. [l] pictorially represents a hadron-coUider event, where a tth final state is produced and evolves by 
including effects of QCD bremsstrahlung in the initial and final state, the underlying event, hadronisation 
and, finally, the decays of unstable hadrons into stable ones. Event generators usually rely on the fac- 
torisation of such events into different well-defined phases, corresponding to different kinematic regimes. 
In the description of each of these phases different approximations are employed. In general the central 
piece of the event simulation is provided by the hard process (the dark red blob in the figure), which 
can be calculated in fixed order perturbation theory in the coupling constants owing to the correspond- 
ingly high scales. This part of the simulation is handled by computations based on matrix elements, 
which are either hard-coded or provided by special programs called parton-level or matrix-element (ME) 
generators. The QCD evolution described by parton showers then connects the hard scale of coloured 
parton creation with the hadronisation scale where the transition to the colourless hadrons occurs. The 
parton showers model multiple QCD bremsstrahlung in an approximation to exact perturbation theory, 
which is accurate to leading logarithmic order. At the hadronisation scale, which is of the order of a 
few Aqcd, QCD partons are transformed into primary hadrons (light green blobs) by applying purely 
phenomenological fragmentation models having typically around ten parameters to be fitted to data. 
The primary hadrons finally are decayed into particles that can be observed in detectors. In most cases 
effective theories or simple symmetry arguments are invoked to describe these decays. Another impor- 
tant feature associated with the decays is QED bremsstrahlung, which is simulated by techniques that 
are accurate at leading logarithmic order and, eventually, supplemented with exact first-order results. A 
particularly difficult scenario arises in hadronic collisions, where remnants of the incoming hadrons may 
experience secondary hard or semi-hard interactions. This underlying event is pictorially represented by 
the purple blob in Fig.[T] Such effects are beyond QCD factorisation theorems and therefore no complete 
first-principles theory is available. Instead, phenomenological models are employed again, with more 
parameters to be adjusted by using comparisons with data. 
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Modern event generators 

The most prominent examples of event generators are the highly successful, well-established programs 
Pythia m and Herwig [2]. They have been constructed over the past decades alongside with exper- 
imental discoveries and most of the features visible in past and present experiments can be described 
by them. However, the need for higher precision to meet the challenges of new energy scales occuring 
at the LHC, the complexity of final states at those scales, the necessity of maintenance and the wish 
to easily implement new physics models have demanded those codes to be rewritten in a modern pro- 
gramming language providing a higher level of modularity. Object-oriented frameworks meet the latter 
requirements, and owing to the community's preference for C++, the new generation of event generators 
is constructed in this programming language. This led to improved re-implementations in form of the 
programs Pythia 8 [S] and HERWIG++ |lj - the successors of the Fortran versions mentioned above - 
and to the construction of the Sherpa event generator jB] . 

In conjunction, in the past decade codes for next-to-leading order calculations have been made available 
to the public; prominent examples include MCFM [6] and NloJET-i-i- [7]. Furthermore, methods have 
been proposed for the consistent matching of next-to-leading order corrections with parton-shower al- 
gorithms [HllSj- Corresponding methods are implemented for example in Mc@NL0 |10j, which is based 
on the Fortran version of HERWIG, in HERWIG [TT] and in some more specialised programs [T^ . 
However, the full next-to-leading order calculations underlying these new techniques are very complex 
and challenging and until today only processes with up to five external legs are under control. On the 
other hand, many important experimental signatures rely on final states with higher multiplicities. This 
has triggered substantial activity in perfecting techniques and tools at tree-level accuracy, such that by 
now several codes are available that can compute corresponding cross sections and generate events in 
a fully automated way. The most prominent examples include Alpgen [T^], CompHEP/CalcHEP jl4j . 
HelaC-PhegaS [H], MadGraPH 16 , WhizaRD [n] and AMEGIC++ [H]. Currently only AMEGIC++ is 
part of, and integrated in, a full-fledged event generator, namely the Sherpa framework. In order to 
translate the multi-particle parton-level events, which are provided by these tools at leading order, into 
hadron-level events, several algorithms have been developed, all aiming at preserving the logarithmic 
accuracy of the parton shower and supplementing it with the exact perturbative leading order result for 
given jet multiplicities. In [TH] an algorithm achieving this goal in e+e^ annihilations into hadrons has 
been presented and it has been extended to hadronic collisions in |20j . A similar algorithm for the dipole 
shower has been discussed in [21], whereas a more different one has been published in [H]. All these 
approaches have been compared in [331IH] and a good agreement has been established. 

1 Sherpa's event generation framework 

Sherpa [S] is an acronym for "Simulation of High Energy Reactions of PArticles" . The program is a 
complete event generation framework that has been constructed from scratch and entirely written in the 
modern, object oriented programming language C++. 

Construction paradigm 

The construction of SHERPA has been pursued in a way largely defined by the following three paradigms: 

• Modularity. Different physics aspects are implemented in almost independent modules, relying on 
a small number of framework and support modules, like, e.g., the event record etc.. Modularity 
allows, for example, to have more than one matrix-element generator or parton shower in parallel, 
with the user be in charge of making a choice. The central module, SHERPA, steers the interplay 
of all other parts and the actual generation procedure. 

• Bottom-to- Top. Physics modules are typically developed in their own right, being tested and 
validated before they are incorporated into the full event generation framework. This in turn 
results in a quite flexible, minimal structure underlying the organisation of event generation. 

• Separation of interface and implementation. In order to facilitate the two requirements above, 
Sherpa relies on a structure where the (nearly independent) physics modules are accessible only 
through physics-specific handlers. These handlers assist SHERPA in generating the event at different 
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stages, each of which is steered through a specific implementation of Event_Phase_Haiidler, such as 
Signal_Process or Jet_Evolution. An example for such an interplay of event phase and physics 
handler is the Matrix_Element_HELndler, enabling the generation of parton-level events either by 
the built-in hard-coded matrix elements or by the matrix-element generator AMEGIC++. This is 
relevant for two event stages, the generation of the signal process and owing to the multijet merging 
procedure the evolution of the jets. 

This overall structure fully reflects the paradigm of Monte Carlo event generation by factorising the 
simulation into well-defined, almost independent phases. Accordingly, each Event_Phase_Handler en- 
capsulates in an abstract way a diff'erent aspect of event generation for high-energy particle reactions. 
This abstraction is then replaced by real physics using handlers, which ensure that the overall event 
generation framework can be blind to the finer details of the underlying physics and its implementation 
in form of a physics module. 

Physics modules 

In the following the main modules currently distributed with Sherpa will be listed and briefly described. 
For a more in-depth discussion, the reader will be referred to either the corresponding section in this 
paper or to the original literature. 

• AmegiC-i-i- [18] . 

This is Sherpa's default matrix-element generator based on Feynman diagrams, which are trans- 
lated to helicity amplitudes using the methods of [51]. It will be described in detail in Sec. [31 
AMEGIC++ has been thoroughly tested for multiparticle production in the Standard Model [2^ . 
Its MSSM implementation has been validated in [27]. AMEGIC++ employs the Monte Carlo phase- 
space integration library PHASIC. For the evaluation of initial-state (laser backscattering, initial- 
state radiation) and final-state integrals, the adaptive multi-channel method of [H] is used by 
default together with a Vegas optimisation J55] of the single channels. In addition, final-state inte- 
gration accomplished by RAMBO [30J and HAAG [SI] is supported. 

• APACIC+-K [321133]. 

APACIC++ generates initial- and final-state parton showering. The shower evolution is governed 
by the DGLAP equations and is ordered in parton virtualities. Coherence effects are accounted 
for by explicit ordering of the opening angles in subsequent branchings. All features needed for 
a consistent merging with matrix elements [191120] are included, however, the main part of the 
merging procedure is implemented in the Sherpa module itself. 

• AMISIC++ [31]. 

This module simulates multiple parton interactions according to [3S] . In SHERPA the treatment of 
multiple interactions has been extended by allowing the simultaneous evolution of an independent 
parton shower in each of the subsequent collisions. This shower evolution is handled by APACIC++. 

• Ahadic. 

AhaDIC is Sherpa's hadronisation package for translating the partons (quarks and gluons) into 
primordial hadrons. The algorithm is based on the cluster-fragmentation ideas presented in [36] , 
which are also implemented in the HerwiG event generators. It should be noted that AhaDIC, 
essentially based on [37], indeed differs from the original versions, cf. Sec. 17.21 

• Hadrons [35] . 

Hadrons is the module for simulating hadron and r-lepton decays. The resulting decay products 
respect full spin correlations (if desired) . Several matrix elements and form- factor models have been 
implemented, such as the Kiihn-Santamari'a model or form-factor parametrisations from Resonance 
Chiral Theory for the r-leptons and form factors from heavy quark effective theory or light-cone 
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sum rules for hadron decays. For further details, see Sec. [51 

• Photons [55] . 

The Photons module holds routines to add QED radiation to hadron and r-lepton decays based 
on the YFS algorithm [40] . The structure of PHOTONS is designed such that the formalism can be 
extended to scattering processes and to a systematic improvement to higher orders in perturbation 
theory, cf. f35] . The application of PHOTONS therefore fully accounts for corrections that are usually 
added by the application of PHOTOS [H] . 

Sherpa itself is the steering module that initialises, controls and evaluates the different phases during 
the process of event generation. Furthermore, all routines for the combination of parton showers and 
matrix elements, which are independent of the specific parton shower are found in this module. For 
details concerning the implementation of these routines, see Sec.[SJ 

In addition to the main modules of Sherpa, there is a set of tools providing basic routines for event 
generation, general methods for the evaluation of helicity amplitudes, some generic matrix elements, etc.. 
Interfaces to commonly used structures like the LhaPDF package [12], or Slha inputs [13] exist as well 
as interfaces to most frequently used output formats like Hepevt and Hepmc. 

In the following the different phases of event generation and the respective realisation within the SHERPA 
program are discussed in detail. Some examples are presented to illustrate the capabilities of the code 
and highlight its special features. 

2 Parton distributions and beam spectra 

Within Sherpa the incoming beam of particles is resolved in two steps. 

• Bcamstrahlung, which leads to interacting photons in the hard process or photon generation 
through laser backscattering is described through appropriate spectra. 

• The parton content of photons and hadrons is parametrised by appropriate parton distributions. 

2.1 Beam spectra 

In the context of leptonic incoming beams, e.g. at a future linear collider, photon beams can be prepared 
through laser backscattering off the (potentially polarised) electrons at very high luminosity 44 . The 
respective scenario can be simulated in SHERPA through a parametrisation that depends on the energy 
and polarisation of the incoming lepton beams as well as on the concrete laser specifications 45 . 
When considering collisions of hadronic beams, e.g. at the Fermilab Tevatron or the CERN LHC, there 
exists a two-photon component of the cross section, which can reasonably be described in the framework 
of the equivalent-photon approximation (EPA) . The EPA relates the hadronic cross section to the inter- 
action cross section of real photons through a two-photon luminosity function 1461 . The corresponding 
implementation in SHERPA has been described in some detail in |47j . 

2.2 Parton distributions 

Sherpa provides a variety of interfaces to standard PDF sets, used to parametrise the parton content 
of incoming particles in hadron collisions. The most commonly employed LHAPDF package [12] is 
interfaced such that the full wealth of PDF sets contained in this package is made available throughout 
the code. 

If photons generated through laser backscattering or radiated off hadrons or nuclei are themselves re- 
solved, their quark and gluon content may be parametrised by a photon PDF, see for example [48] . In 
this context the event generation is essentially equivalent to the procedure for hadronic initial states, 
although with varying centre-of-mass energy. This can easily be understood, since on the perturbative 
side of the simulation it only matters that there is a PDF. No information about the incoming beam 
being a photon is needed. 
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Fig. 3 



Translation of a Feynman diagram into a 
phase-space parametrisation. Ds^a denote 
symmetric or asymmetric decays - the lat- 
ter ones reproduce the typical feature of 
coUinear emissions of particles notorious 
for gauge theories with massless spin-one 
bosons. The propagator terms for massless 
particles, Po, peak at the minimally allowed 
invariant mass. 



Fig. 



2 Generic form of a three-point interaction 
vertex. The interaction of particles 1, 2 and 
3, with 1 incoming while 2 and 3 are out- 
going, is defined by the left- (cl) and right- 
handed (c_r) coupling, the SU{3) colour op- 
erator T, and the Lorentz structure of the 
interaction L. 



3 Hard matrix elements and phase-space integration 

Matrix elements for the hard scattering process can be provided either by a fast internal library of 
hard-coded 2^2 processes or in an automated way by using the matrix element generator AMEGIC++. 
AMEGIC++ [IH], acronym for "A Matrix Element Generator in C++", is a multi-purpose parton-level 
generator. It is a convenient tool for calculating cross sections of nearly arbitrary scattering processes 
at tree level based on e.g. the following models: 

• Standard Model, 

• Extension of the SM by a general set of anomalous triple and quartic gauge couplings [49] , 

• Extension of the SM by a single complex scalar [50], 

• Extension of the SM by a fourth generation, 

• Extension of the SM by an axigluon |51j . 

• Two-Higgs-Doublet Model, 

• Minimal Supersymmetric Standard Model, 

• ADD model of large extra dimensions [55] . 

Besides calculating production and decay rates AMEGIC++ is used to generate parton-level events vifithin 
the event simulation framework of Sherpa. Therefore, these partonic events can easily be supplemented 
by parton showers and linked to the hadronisation, which in turn yields realistic hadronic final states. 

Given a physics model and a certain process AMEGIC++ automatically generates the corresponding set 
of tree-level Feynman diagrams from the complete set of interaction vertices possessed by the model0 
A vertex is thereby defined by the incoming and outgoing particles, the left- and right-handed coupling, 
the explicit SU{S) colour structure of the interaction, and the related Lorentz structure. For a schematic 
representation of a generic three-point function, see Fig. [2l Four-point interactions have one more 
outgoing particle, otherwise are set up completely equivalently. 

The available colour operators T are 1, Sij , Sab, T^j , fate and products of those, e.g. fabefcde- Some vertices 
are decomposed into several colour or Lorentz structures, for example the 4-gluon vertex is represented 
by three products of colour structures with a respective Lorentz operator. Because of this quite general 
and explicit treatment of colour new physics models can be implemented in a straightforward manner 
even though they might imply interactions with colour structures absent in the SM. 
The Feynman diagrams constructed from the vertices are translated into helicity amplitudes relying on 
a formalism similar to the one described in [25] and extended to include also spin- two particles in [53] . 
During this translation process the chained Lorentz structure of each single diagram is mapped onto 

^ Note that unitary gauge is considered for the vertices; moreover AMEGIC++ is currently limited to three- and four-point 
interactions. 
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suitable products and sums of hclicity-amplitude building blocks. The diagrams are then grouped into 
sets of amplitudes with a common colour structure. Based on them, the exact matrix of colour factors 
between amplitudes is explicitly calculated using SU (3) algebra relations. A number of refinements 
of the hclicity method have been implemented in the code to eventually speed up the matrix-element 
evaluation and cope with issues that arise when dealing with extensions of the SM. Concerning the latter, 
cf. Sec. KTT\ 

For external massive or massless gauge bosons, explicit polarisations are enabled. This allows to calculate 
polarised cross sections. Furthermore, the numerators of spin-one and spin-two propagators can be 
replaced by sums over suitably defined polarisations for off-shell particles, thereby disentangling nested 
Lorentz structures emerging for amplitudes with many internal bosons. As a result, the generic 3- and 
4-point Lorentz structures of the considered physics model are the only basic helicity-amplitude building 
blocks, which are needed by AMEGIC++ to construct arbitrary processes. The complete sets of vertex 
structures appearing in the SM, the MSSM and the ADD model, plus a comprehensive list of generalised 
triple and quartic gauge-boson interactions have been implemented. Fortunately, these sets also form the 
basis for many other extensions of the Standard Model and, hence, such BSM models can be incorporated 
without the need to provide new helicity-amplitude building blocks. 

Employing the helicity formalism allows to accelerate the matrix-element evaluation by making use of 
symmetries among different Feynman graphs. One example is given by diagrams with equal colour struc- 
ture that have common factors, as it is illustrated in Fig. HI Accordingly, identical sub-amplitudes can be 
factored out and the cut graphs can be added, thereby reducing the number of complex multiplications 
to be carried out. The optimised helicity amplitudes for the particular process are stored in library files. 

The second main task, which needs to be accomplished by AMEGIC++, is the generation of suitable phase- 
space integrators, which allow an efficient cross-section evaluation and subsequent event generation. The 
employed algorithm combines process-specific a priori knowledge about the integrand with self-adaptive 
Monte Carlo integration techniques [28l[29l[54] ■ 

The Feynman diagrams for a given process are analysed to create specific phase-space mappings, which 
generate non-uniform momentum distributions with weights that approximate the corresponding (squared) 
diagram. The phase space for a single diagram is parametrised using invariant masses of the diagram's 
propagators and variables related to the angles in particle splittings. As exemplified in Fig.[2]the mapping 
(and the corresponding weight function) can be assembled out of a few generic building blocks generating 
corresponding weighted distributions of those parameters. As the matrix elements themselves, all phase- 
space mappings are stored in library files. For the integration of a scattering process, all contributing 
channels are combined in a self-adaptive multi-channel integrator, which automatically adjusts to the 
relative importance of the single phase-space maps to minimise the variance. 

The efficiency of the integrator is further improved by applying the self-adaptive VEGAS algorithm on 
single phase-space maps. VEGAS [29J is very efficient in the numerical adaptation to functions, which 
are factorisable into a product of one-dimensional functions. Although this is clearly unlikely for a full 
matrix element, the structure represented by a single phase-space mapping fulfills this condition. VEGAS 
thus allows an adaptation to more detailed structures generated e.g. by phase-space cuts or the spin of 
the particles. These structures clearly are beyond the approximations inherent in the mappings. 
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Both the hcUcity ampUtudcs and the phase-space parametrisations, are semi-automaticall50 compiled 
and hnked to the code before the actual integration and event generation can start. 
In recent versions of Sherpa several improvements to the algorithm outlined above have been imple- 
mented. They shall briefly be discussed here. 

• Besides calculating the full set of Feynman diagrams for a given final state it is now possible to 
evaluate only certain resonant graphs from decays of unstable intermediate particles. Accordingly, 
the process is decomposed into the actual production and the decay (s) of the resonance(s). For a 
number of intermediate states, labelled i, the combined amplitude reads 

^(") = <,:r^« n ^^-^t!- (1) 

i (E decays 

The (E) symbol represents the colour and spin correlations between the production and decay am- 
plitudes. The factor Vi accounts for the intermediate propagator. Two options are realised: either 
the resonances are forced on their mass shell, or the full off-shell propagator is considered. Both 
approaches fully preserve colour and spin correlations, however, only the first method yields strictly 
gauge-invariant results, owing to the usage of on-shell matrix elements for production and decays. 
Although gauge invariance cannot be guaranteed when the full off-shell propagator is considered, 
this ansatz naturally incorporates finite-width effects. 

This decay treatment can be applied iteratively, therefore, provides a very elegant way to simulate 
entire decay chains as they appear for example in supersymmetric theories. Furthermore, the 
treatment allows for an easy incorporation of rt-body decays as the decay amplitudes ^doc are not 
restricted to 1 — > 2 processes. 

• Despite the optimisation strategies employed in 

AMEGIC++ the evaluation of Feynman diagrams for purely strong-interacting processes, in par- 
ticular n-gluon scattering, is still quite cumbersome. For this case, the matrix elements can be 
evaluated using Cachazo-Svrcek-Witten vertex rules [55j, which are based on maximal helicity 
violating (MHV) amplitudes as building blocks. For small and intermediate particle multiplicities, 
typically up to 2 — > 5, these rules give the most compact expressions, cf. |56) . 

Algorithmically there is no upper limit on particle multiplicities for which matrix elements can be gener- 
ated. Practically, calculations are however limited by the accessible computer power. Another important 
criterion is given by comparing Amegic++'s performance to alternative approaches for matrix- element 
generation, e.g. [13 1 ? J57j . Using the Feynman diagrammatic approach (including the phase-space integra- 
tion method described above), AMEGIC++ performs reasonably and can compete for scattering processes 
that typically involve up to a few thousand diagrams. In the framework of the Standard Model this is 
sufficient for processes with a total number of up to eight or nine particles, or a maximum of six partons 
in the context of evaluating hadronic cross sections. The latter limit is not only due to the proliferation 
in the number of Feynman diagrams but also the number of different parton-level processes contributing 
to a jet cross section, which typically are all desired to be taken into account at once. Accessible particle 
multiplicities can of course be increased by employing the above mentioned decay treatment, since it 
restricts the number of amplitudes by accounting for certain resonant graphs only. 

The built-in matrix-element generator AMEGIC++ renders Sherpa a very powerful tool to study hard 
production processes, in particular many-particle final states. The underlying tree-level calculations 
thereby naturally account for aspects of multi-particle production processes that often are only approx- 
imated in alternative approaches. Thus, a clear strength of complete matrix-element calculations is 
the proper incorporation of off-shell (finite-width) effects and quantum interferences between different 
diagrams contributing to the same final state. The correct treatment of angular correlations between 
final-state particles can be achieved as well. All of these aspects are important in approaching a realistic 
description of both signal and background processes. 

■^In detail this means that after a first run generating all process libraries, these libraries have to be compiled and 
dynamically linked to the main code by executing the automatically generated "makelibs" script. 
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Number of jets 


e+Ue. + QCD jets 





1 


2 
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AMEGIC++ 
Comix 


5432(5) 
5434(5) 


1279(2) 
1274(2) 


466(2) 
465(1) 


185.2(5) 
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77.3(4) 
77.5(3) 




a [pb] 


Number of jets 


e"e+ + QCD jets 
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2 


3 


4 


AMEGIC++ 
Comix 


723.0(8) 
723.5(4) 


188.2(3) 
187.9(3) 


69.6(2) 
69.7(2) 


27.21(6) 
27.14(7) 


11.1(1) 
11.09(4) 





Number of jets 
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2 


3 


4 
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AMEGIC++ 
Comix 


331.0(4) 
331.0(4) 


22.78(6) 
22.72(6) 


4.98(1) 
4.95(2) 


1.238(4) 
1.232(4) 



a [pb] 


Number of jets 


tt + QCD jets 





1 


2 


AMEGIC++ 


754.4(3) 


747(1) 


520(1) 


Comix 


754.8(8) 


745(1) 


518(1) 



Tab. 1 Comparison of Standard Model production cross sections at the LHC using the MC4LHC parameter 
setup [5T]. In parentheses the statistical error is stated in units of the last digit of the cross section. 



3.1 Standard Model production processes 

All generic interactions stemming from the Standard Model Lagrangian have been implemented, allowing 
the generation of matrix elements for arbitrary processes involving all three generations of quarks and 
leptons, all SM gauge bosons and the Higgs boson. 

In addition to the generic SM interactions effective interaction vertices between a Higgs boson and 
massless gauge bosons are available in AMEGIC++. The coupling to gluons is mediated by a top-quark 
loop and modelled through the effective Lagrangian [55] 

^'ggH 9ggH-^;^G^^G'^'^H , (2) 

where G^^ = S^A" — — gf°-'"^A''^A'^ is the gluon field-strength tensor. The effective coupling gggH 

can be calculated either for finite top-quark mass or in the limit mt oo. 

For the coupling of a Higgs boson to photons, mediated by top-quark and W-boson loops, the effective 
Lagrangian reads 

^ffH = ^F^^F'^'-H, (3) 
where F^j^i, = 9^/1° — dtyA"^ is the photon field-strength tensor. 

The validation of Sherpa's core - providing the matrix-element calculations - is an important aspect in 
the development of the Sherpa Monte Carlo generator. AMEGIC++ has successfully been tested against 
various other programs for a great variety of processes. This includes cross sections and distributions for 
photon-photon collisions [59], e"'"e~ annihilations at different centre-of-mass energies [60"26', and particle 
production in hadron-hadron collisions ^Ql^. In Tab. [Tj exemplary results of a cross-section comparison 
for some LHC key processes are presented. For this comparison between AMEGIC++ and the independent 
matrix-element generator COMIX, the calculational setup used in pT] has been employed. The results of 
both codes agree very well within the statistical uncertainties, indicated by the numbers in parentheses, 
proving the correctness of the respective calculations. 

To illustrate the relevance of accounting for spin correlations we consider the process 77 W^W~ 
IV +]fiT at the CERN LHC, either in the continuum or mediated by a Standard Model Higgs boson 
of mass ran — 160 GeV. Figure shows the azimuthal separation of all possible lepton combinations 
emerging from the VF-decays. Apparently, this observable provides an excellent handle on the irreducible 
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Fig. 5 Lepton azimuthal separation in 77 W^W II' +fJT with and without an intermediate Standard 
Model Higgs boson of mn = 160 GeV. 



W~^W~ background process. Note that both spectra have been normahsed to unity in order to compare 
their respective shapes. To parametrise the initial-state photon energy, the spectra described in Sec. [D 
have been employed. 

3.2 Matrix elements for physics beyond the Standard Model 

The implementation of new physics models into AMEGIC++ proceeds in two steps. First, the new particles 
and the parameters of the model have to be declared and corresponding setting routines need to be 
providedlfl Second, the generic interaction vertices of the theory have to be defined. For each new three- 
and four-point function, this includes the specification of the incoming and outgoing particles and the 
couplings, as well as the assignment of the SU (3) colour structure of the interaction and its Lorentz 
structure. Both steps are centralised in the module MODEL such that the parametric quantities of the 
model can consistently be used in all Sherpa modules. 

The only limitation of the AMEGIC++ approach, concerning its extension to new physics scenarios, is the 
number of Lorentz structures or helicity-amplitude building blocks known to the program. Behind each 
Lorentz structure there resides a corresponding calculational method for the evaluation of a respective 
sub-amplitude. Although implementing a new building block is a well defined task, users that wish to 
do so should contact the authors for advice and help. However, there already exists a large variety of 
interaction operators for physics beyond the SM. In the following the new physics models available in 
AMEGIC++ are briefly reviewed and some implementation details will be highlighted. 

3.2.1 The MSSM implementation 

For the Feynman rules of the i?-parity conserving Minimal Supersymmetric Standard Model, the con- 
ventions of [62j are used. The general set of interaction vertices derived there, and as such implemented 
in AMEGIC++, include not only a full inter-generational mixing of the squark and slepton fields but also 
permit the inclusion of CP violating parameters. Furthermore, they include finite masses and Yukawa 
couplings for all three fermion generations. 

Necessary ingredients when dealing with the MSSM are specific Feynman rules for Majorana fermions or 
fermion number violating interactions. To unambiguously fix the relative signs among Feynman diagrams 
involving Majorana spinors the algorithm described in [63j is used and the general set of fermion Feynman 
rules given therein has been implemented. Accordingly, the explicit occurrence of charge-conjugation 
matrices is avoided, instead a generalised fermion flow is employed that assigns an orientation to complete 

^ Essentially this is accomplished by reading numerical values of masses and other parameters from input files. 
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cr(W+W- ^ 


X) [fb] 






Final 


MadGraph 


Whizard 


Sherpa 


state 


y/s =0.5 TeV 


y/s =2 TeV 


y/s =0.5 TeV 


y/s =2 TeV 


y/s =0.5 TeV 


y/s =2 TeV 


xUI 


3.8822(2) 


1.2741(4) 


3.8824(1) 


1.27423(8) 


3.8821(2) 


1.2741(1) 


x?x^ 


121.29(1) 


24.47(1) 


121.2925(7) 


24.472(3) 


121.296(6) 


24.477(1) 


x?x^ 


6.8936(7) 


12.880(7) 


6.8934(2) 


12.8790(8) 


6.8938(3) 


12.8793(6) 


x?x^ 


1.4974(1) 


9.707(5) 


1.4973(6) 


9.7064(7) 


1.49735(7) 


9.7078(4) 


x^x^ 


5996.5(4) 


1.0415(6)e3 


5996.57(2) 


1.04150(5)e3 


5996.4(3) 


1.04148(5)e3 


x^x^ 


— 


365.6(2) 


— 


365.615(6) 


— 


365.63(2) 


x^x^ 


— 


467.8(2) 


— 


467.775(8) 


— 


467.77(2) 


nn 


— 


82.35(3) 


— 


82.347(8) 


— 


82.352(4) 




— 


138.20(5) 


— 


138.18(1) 


— 


138.205(7) 


xlxl 




117.78(4) 




117.80(1) 




117.786(6) 


xtxi 


3772(1) 


944.3(8) 


3771.6(4) 


944.2(1) 


3771.8(2) 


944.32(5) 


xtx2 




258.3(2) 




258.37(4) 




258.36(1) 


xtx2 




131.0(1) 




130.98(2) 




130.966(7) 



Tab. 2 Sample cross sections (in fb) calculated using MadGRAPH, Whizard and SherpA for the production 
of neutral and charged gauginos in WW scattering at fixed centre-of-mass energies of 0.5 and 2.0 
TeV. The numbers in brackets reflect the absolute statistical error on the last digit. The consid- 
ered supersymmetric spectrum corresponds to the SPSla benchmark scenario [Bl], generated with 
SOFTSuSY [65]. Additional details on the calculational setup can be found in [27] . 



fermion chains. This uniquely determines the external spinors, fermion propagators and interaction 
vertices involving fermions. Furthermore, negative mass eigenvalues of the physical neutralino fields 
are taken into account at face value in the propagators and spinor products in the helicity-amplitude 
expressions. In this way a redefinition of the neutralino fields and couplings can be avoided. 
Files that are conform with the SUSY Les Houches Accord (SLHA) [l^ are used to input the MSSM weak- 
scale parameters. Such files are provided by external spectrum calculators, see for instance Ref. [65]. The 
parameters are translated to the conventions of |62j and accordingly fix the couplings of all interaction 
vertices. Note that so far Sherpa only supports SLHA-vl and the generalisations presented in [HS] have 
not been implemented yet. This means that further assumptions on the physical model are implicit. The 
mixing of the squark and slepton fields is restricted to the third generation only and all SUSY breaking 
parameters and the Yukawa couplings are assumed to be real such that the model realises exact CP 
symmetry. 

To verify the correctness of the implemented Feynman rules and the translation of the SLHA inputs to 
the conventions of [^2] , detailed comparisons against the matrix-element generators MADGraPH [TB] and 
Whizard [T7] , have been carried out [27] . Several hundred cross sections for supersymmetric processes 
have been compared and mutual agreement between the three independent codes has been found. As 
an example some of these results are listed in Tab. [5] Besides comparing explicit cross sections, thereby 
testing all phenomenologically relevant couplings, various unitarity and gauge-invariance checks were 
performed and supersymmetric Ward- and Slavnov- Taylor identities were tested. 

With the complete MSSM Lagrangian being available in a full-fledged matrix element generator like 
AMEGIC++ detailed aspects of SUSY phenomenology can be studied. Examples include the impact of 
off-shell kinematics on MSSM signals [67]l68j . interference effects between various SUSY processes leading 
to the same final state (irreducible backgrounds) 23 > oi' the impact of additional hard QCD radiation 
[691170] . Furthermore, non-trivial production processes such as the weak-boson-fusion production of 
supersymmetric particles can easily be considered [71] . 

3.2.2 The ADD model of large extra dimensions 

The Arkani-Hamed-Dimopoulos-Dvali model [52J extends the usual 3-|-l-dimensional SM by 5 additional 
compactificd spatial dimensions where only gravity can propagate. In this model the usual 4-dimensional 
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Fig. 6 Dependence of the total cross section on the ADD scale Ads for jet-graviton production. The SM 
background at the LHC is also shown. A phase-space cut of p™"''' > 1 TeV is imposed. The dashed 
lines are for Mq^ < Ms, while the solid lines are for y/§ < Ms- 



Planck scale is related to a fundamental scale Ms by 

M|i = SnR^Mp-^ (4) 

where is the volume of the compactified extra-dimensional space. A sufhciently large compactification 
radius R allows Ms of the order of 1 TeV, providing an explanation of the mass hierarchy between 
the electroweak and the (4-dimensional) Planck scale. This has observable consequences for TeV-scale 
colliders. An advantage of this model is that it introduces only two new parameters, the scale Ms and 
the number of extra dimensions 6. 

In a 4-dimensional picture the ADD model leads to the coupling of an abundance of Kaluza-Klcin 
graviton states to SM fields. Feynman rules for those interactions can be derived from a linearised 
gravity Lagragian [721173], 

The AMEGIC++ implementation includes all 3- and 4-point interactions between gravitons and SM fields. 
It features the generation of matrix elements, which include the exchange of virtual gravitons and the 
production of real gravitons. To obtain physical cross sections this model additionally requires a sum- 
mation (usually replaced by an integration) over the Kaluza-Klein states of the graviton. For the two 
cases this requires different approaches: 

• Virtual graviton production: 

The integration is performed analytically at the level of the graviton propagator. A number of 
parametisations appeared in the literature, which include a necessary cut-off for masses above 
some scale of the order of M5. Three widely used conventions have been implemented [T ^ FfS l IT^. 

• Real graviton emission: 

The sum over accessible Kaluza-Klein states is evaluated together with the phase-space integral 
using Monte Carlo techniques. 

For details of the implementation, cf. [53] . 

A simple example for the calculation of an ADD cross section with AMEGIC++ is given in Fig. which 
presents the LHC monojet cross section from ADD graviton production for various numbers of extra 
dimensions as a function of the scale Ms in the notation of [73]. Figure [6] also shows the SM background 
from Z bosons decaying into a neutrino pair0 

All new Lorentz structures implemented for the ADD model are completely generic. This allows to use 
them in other models that feature couplings to gravitons. 

^Note that the results presented in Fig.[6]can be compared to those of Fig. 5(a) in |75| . 
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3.2.3 Anomalous interactions of the SM gauge bosons 

Besides the complete set of SM weak-gauge boson interactions AMEGIC++ includes a number of effective 
Lagrangians describing anomalous triple and quartic gauge interactions. To account for the missing UV 
completion of the effective theory approach, a simple unitarisation method can be applied. 

• WWV interactions: 



The general set of operators describing the interaction of two charged vector bosons and a neutral 
one, 

Cwwv/gwwv = igYiWl^W^V -WlV.Wn 



m 



w 



-glWlW^id^^V + d'^V^) + gXe^'^P'^iWl dpW,)V„ 

+ V + ^^"''^WI^W::V,„ , (5) 

cf. [76], is implemented. Here V^^ denotes either the photon or the Z field, is the W~ field, 
W^, = d^W, - d,W^, V^, = d^V, - d,V^, = \e^upaVP" and {aXb) = A{d,,B) - {d^,A)B. 
The overall coupling constants are given by 

51^1^7 = — e and gwwz = —ecotOw ■ (6) 

The SM values of the individual couplings are gf^'^ = kz/-^ = 1 while all others vanish. 

• Quadruple interactions: 

The following SU (2) custodial symmetry conserving quartic interactions are available: 

A = a^e^ (iwlW^^'W.W + ^iWlW^^f + -^WlZ^W.Z" + -^{Z''Z'')A , (7) 

A = aJ{WlWPf + ^WlZ^'W,Z'' + -^{Z^'Z^'f\, (8) 
cf. |77j . In the SM limit the parameters and as are identical zero. 

• j-Z interactions: 

The general anomalous coupling of two on-shell neutral gauge bosons (Vi and V2) to an off-shell 
boson (V3) has been given in ^6\. The corresponding Feynman rule is displayed in Fig. [T] The 
vertex functions FviVaVg* for V2 — Z and V2 = 7 are given in Eqs. (O and ([TU]) . respectively. In the 
SM all coupling parameters, i.e. fY, fY 7 ^ equal zero. 



^ Z-yVyi^i<l2,P) - 2 \'^l\l2g -923 j 

-P" [{Pq2)g^^ - q^^P^] + 4e''"''''92p + ^P" e^'^'" Ppq2a\ 
: ^z ) 

It should be noted that the most general anomalous coupling between three off-shell neutral gauge 
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Via{qi 




^2/3(92) 

Fig. 7 Feynman rule for the anomalous triple neutral gauge boson vertex. 



bosons allows more coupling terms [78], which are, however, not implemented in the current ver- 
sion. To account for this a symmetrised version of the above vertex is used outside the on-shell limit 
of two of the vector bosons. For example, the interaction of three off-shell Z bosons is modelled 

through Tz'Z'Z* = ^ZiZ2Z* + T^aZaZ' + ^ZsZiZ*- 



Unitarisation: 

Owing to the effective nature of the anomalous couplings unitarity might be violated for coupling 
parameters other than the SM values. For very large momentum transfers, such as those probed at 
the LHC, this will lead to unphysical results. As discussed in [7S] this can be avoided by introducing 
energy dependent form factors for the deviation of coupling parameters from their Standard Model 
values: 

= 7-^ ^° . „, (11) 

^ ' (1 + s/A2)» ^ ' 

where s is the partonic centre-of-mass energy, A represents the new physics scale, and the exponent 
n is typically chosen to be n = 2. 



3.2.4 Further BSM models 

Besides the major model implementations described above AmegiC-i-i- offers various other physics sce- 
narios that shall briefly be listed here: 

• Two-Higgs-Doublet model: 

As byproduct of the MSSM implementation the THDM of type II can be studied. In contrast to 
the MSSM, the relevant input parameters of the extended Higgs sector (masses of the Higgs bosons 
h^, H^, A^, and , the mixing angle between the two doublets, a, and the ratio of their vevs, 
tan /?) have to be specified explicitly. 



• Phantom-Higgs model: 

This model has been discussed in [5D]. It emerges by adding a complex scalar gauge singlet to 
the SM, which interacts with the SM particle through a mixed quartic coupling with the original 
Higgs doublet of the SM. After symmetry breaking, a massless pseudoscalar and a massive scalar 
augment the original SM particle spectrum. The latter will undergo a mixing with the original 
Higgs boson, therefore dilute its couplings to the other SM particles through a mixing parameter. 
On the other hand, the massless pseudoscalar does not interact with the SM world, apart from the 
Higgs boson, leading to a sizable invisible partial width of both scalars. The relevant parameters 
in this model are the masses of the two scalars, their mixing angle, tan 6, and the ratio of vevs, tan /3. 



• 4th generation: 

In this model a fourth generation family has been included, i.e. a lepton £4, a neutrino 1^4, an 
up-type and a down- type quark U4 and ^4, respectively. They are parametrised by their masses, 
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which also serve as input parameters. In addition, both may decay to the Standard Model par- 
ticles through the charged current weak interaction where the mixing to the other generations is 
parametrised by one mixing angle each. In the lepton case, this implies that all £4 and 1/4 must 
decay into a r-neutrino or a r-lepton, respectively, while the pattern in the quark case is a bit more 
involved. There the mixing matrix emerges by multiplying the original 3x3 CKM matrix with a 
rotation matrix in the 34-sector, cf. j80j . 



• Axigluons: 

Here, an axigluon, i.e. a massive axial vector with the colour quantum numbers and interactions of 
the gluon is added to the Standard Model, cf. [51]. The only input parameter is the mass of this 
particle. The implications of such a particle have been discussed, e.g., in [§1\ . 

There exist further (though private) implementations of alternative physics scenarios that are not publicly 
available yet but can be obtained from the respective authors on request, see e.g. Refs. |HHI13]. Note 
that to further alleviate the implementation of new physics ideas and to allow for an easy distribution 
an interface to the FeynRules package [84J is currently being developed. 



4 Initial- and final-state radiation 

QCD parton evolution and the occurrence of jets can be understood theoretically when the structure of 
perturbative amplitudes is examined in the kinematical regime of intrajet evolution, i.e. where two or 
more partons get close to each other in phase space. Whenever this happens, any QCD matrix element 
squared factorises into a matrix element squared containing the combined "mother" parton and a univer- 
sal function describing the splitting into the "daughters" . In this limit, the theory becomes semi-classical 
and can be understood in a Markovian approach, where a single initiating parton develops a cascade of 
independent branchings. This is the basic concept of any shower Monte Carlo. Potential differences then 
arise in the factorisation scheme only. The default shower generator within Sherpa is APACIC++ [32.33]. 
It is essentially based on virtuality ordered DGLAP parton evolution with superimposed angluar ordering 
constraints. 



4.1 The parton cascade APAClC-K-h 

Final-state showering in APACIC++ proceeds along the lines of [SS]. evolution of QCD partons is 
simulated by 1 ^ 2 splittings, which occur with differential probability 

d^bfalch.g d / /■*' dt a^M£j)) ^ 1 



Here a denotes the flavour of the splitting parton and Pab{z) are the unregularised Altarelli-Parisi kernels 
in four dimensions for the splitting a —^ be. The transverse momentum k^ is given with respect to the 
axis defined by the direction of the decayer. The condition kj^ > k_L_i„in, with kx ^iin some cutoff value 
for k_L, yields the boundaries of the z- integral, z^in and i„iax, cf. |33j . Evolution and splitting variable 
are defined by 

t =pl-ml , and ^ = ^ ' (13) 

respectively, rua being the on-shell mass of parton a. The splitting variable z is related to the light-cone 
momentum fraction z ^p^/Pa (with the "+" direction defined by pa) through 



2SnE^ + z{ta~ K^El) ^ta + h-t^ where ^=1 + ^1--^. (14) 

Although defined in an apparently non-covariant way, the splitting variable actually is Lorentz invari- 
ant |85j . Colour coherence during evolution is taken into account by an explicit angular veto, which 
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means that a branching is rejected if the opening angle of the emission is larger than the one of the 
previous branching. 

Initial-state showering proceeds in the backward-evolution picture along the lines of [SH] . The differential 
branching probability reads 

where the ratio fb{x/z,t)/ fa{x,t) accounts for the change of parton distributions in each shower step. 
The splitting variable z can be reinterpreted as 

z= i, (16) 
s 

s and s' being the partonic center of mass energies before and after the branching, respectively. This 
immediately yields the relation x' — x/z, thus is partially defining the kinematics after the branching. 
An important issue for DGLAP shower algorithms is the convention to implement kinematic constraints 
once a splitting generates recoil owing to the branching parton going off mass-shell. The recoil strategy 
seems ambiguous because the branching equations are independent of it. In fact, for the derivation of 
the DGLAP equation, it may be assumed that there is a spectator parton aligned along the same axis 
as the splitter, but with opposite direction. This leads to the following approach for APACIC++: 

• In final-state branchings the parton, which originates from the same splitting as the branching 
parton, takes the recoil. If this parton has been already decayed, the decay products are boosted 
accordingly. This amounts to redefining the splitting variable of the respective branching by 



z 



where t and t' denote original and reassigned virtualities, respectively. 

• In initial-state branchings all remaining partons take the recoil. For any splitting b ^ a, the process 
is redefined with parton b rather than a aligned along the beam axis and s s/z. 

The AltarcUi-Parisi splitting functions, Patiz), are taken in the quasi-coUinear limit and include mass 
effects of the emitting parton, cf. [57] . 



Pgg{z) = Cf 



1 



2^i, 



1 - 

Pg,{z) ^Tn[z^ + {l-zf~p.%] , (18) 
z 1 — z 



Pasi^) =Ca 



+ z{l~z) 



1 - z 

Quark-mass dependencies are given in terms of the dimensionless variable 
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{p.+Pj?-^l- ^^^^ 



APACIC++ has been thoroughly tested and validated [33]. Modifications to the shower algorithm, which 
were necessary for a proper implementation of the CKKW merging procedure, have been implemented 
in full generality, see Sec. El 



4.2 Showering off heavy resonances 

APACIC++ is also equipped with the possibility to generate radiation off intermediate heavy resonances, 
once these are described by separable production and decay processes, cf. Sec. [31 An example is top-quark 
pair production, which plays a significant role, both as a signal for a better measurement of Standard 
Model parameters and as a background to new physics searches. 
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The QCD radiation pattern in heavy-flavour decays has been thoroughly investigated in Refs. [55] , 
Within APACIC++ a rather simple strategy is employed. QCD radiation off the decaying heavy particle 
is described by the standard parton shower with massive splitting functions, except for two modifications: 

1. In ordinary final-state parton showering, the mother particle goes off-shell, while the daughters 
retain their respective on-shell masses. In showering off decaying heavy particles, on the contrary, 
the mother particle retains its on-shcU mass, while the daughter of the same flavour goes off-shell 
with decreased virtuality. 

2. The maximally allowed phase-space volume in a branching process of a decaying particle is reduced 
by the factor 



""^^ d$(2)(t,4,i,) X{t,h,t,) ^^^^ 

where A(a, 6, c) = {a — b — c)^ — Abe, t and t' are the virtualities of the decaying particle before 
and after the emission, respectively, and tt and tc are the virtualities of the decay products. This 
correction weight corresponds to a decrease in phase-space volume d$*^^^ owing to a decrease in 
three-momentum |p™j,| in the centre-of-mass frame of the decayed parton. 

To correctly describe the decay process of the heavy flavour, it is vital to respect spin correlations 
between production and decay amplitudes. Within Sherpa, this is done by firstly computing the full 
matrix element for production and decay of the heavy flavour and subsequently adding in the parton 
evolution of the intermediate state. In this respect, it has to be defined, how the kinematics of the 
decay products is to be reconstructed, once a parton emission has occurred in the shower. The following 
strategy is used: 

1. If the decaying heavy particle keeps its mass, i.e. if the radiation occurs in the production part 
of the process, the decay products are simply boosted into the new centre-of-mass frame of the 
decayer. 

2. If the decaying heavy particle does not keep its mass, i.e. if the radiation occurs in the decay part 
of the process, the decay products are reconstructed such that in the centre-of-mass frame of the 
daughter the momenta point into the same direction as they did in the centre-of-mass frame of the 
mother before. 

This procedure largely retains the correlations between the final-state particles initially described by the 
matrix element. 



5 Combining LO matrix elements and parton showers 

One of the most important challenges for simulating events at modern high-energy collider experiments 
is the accurate theoretical description of multijet final states. They constitute the testbed for many 
new-physics searches, whose success will largely depend on our understanding of the Standard Model 
multijet production mechanisms. The best theoretical tools available are therefore needed to approach 
this problem as accurately as possible and obtain reliable estimates of the Standard Model backgrounds 
to experimental analyses. Underlying hard processes should be accounted for by full matrix- element 
calculations and the subsequent evolution and conversion of hard partons into hadronic jets should be 
modelled by QCD parton cascades and phenomcnological hadronisation models. However, there are 
several scales involved that determine the thorough development of an event which makes it difficult to 
unambiguously disentangle the components belonging to the hard process and the hard-parton evolution. 
Given an n jet event of well separated partons, its jet structure is retained when emitting a further 
collinear or soft parton only. An additional hard, large-angle emission, however, gives rise to an extra jet 
changing the n to an n + 1 jet final state. Accordingly a merging scheme for matrix element calculations 
and parton showers has to be defined, that determines on an event-by-event basis which possibility has 
to be followed. The primary goal is to avoid double counting by preventing events to appear twice, i.e. 
once for each possibility, as well as dead regions by generating each configuration only once and using 
the appropriate path. 
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Basically two different solutions for this problem have been proposed and implemented in the past few 
years, the MLM scheme [32], relying on a geometric analysis of the unconstrained radiation pattern 
in terms of cone jets, and the CKKW scheme or CKKW merging algorithm [IHlIinij employing an 
analytical reweighting of the matrix elements supplemented by a constrained parton-shower evolution. 
A reformulation of CKKW to a merging procedure in conjunction with a dipole shower (CKKW-L) has 
been presented in [3T]. Common to all schemes is that sequences of tree-level multileg matrix elements 
with increasing final-state multiplicity are merged with parton showers to yield a fully inclusive sample 
with no double counting and correct at leading-logarithmic accuracy. On the theoretical side there have 
been various studies to compare the different approaches, see e.g. Refs. [5^23p24) . and the first dedicated 
experimental analyses on comparing these improved Monte Carlo predictions with actual data have been 
presented [51111911191^ . However, there is more work needed, on the one hand to understand the systematic 
uncertainties of the various schemes, and, on the other hand further experimental inputs are needed to 
validate the available tools. This is crucial in prospect of the LHC providing first measurements soon. 
The Sherpa generator provides a general and largely process- independent implementation of the CKKW 
algorithm [S^, which has extensively been validated for e~^e~ collisions into hadrons ^93|l33j. the pro- 
duction of single vector bosons at the Fermilab Tevatron [M] and the CERN LHC and for W^W^ 
production at hadron colliders, see [96] . The combined treatment of matrix elements and parton show- 
ers is supplemented by a multiple-interactions description, which respects the jet-production scales of 
the primary process. In the following, these key features of Sherpa shall be discussed in more detail. 
The actual merging algorithm employed will be presented in Sec. 15.11 Aspects of its application will be 
discussed in Sec. 15.21 

5.1 The algorithm implemented in SHERPA 

The idea underlying CKKW is to divide the phase space of partonic emissions according to a kx measure 
[97] into a regime of jet production, described by appropriate matrix elements, and a regime of jet 
evolution, described by parton showering. The separation is defined by the merging scale, denoted by 
Qcut- The matrix elements are reweighted by as coupling factors and terms that arise from analytic 
Sudakov form factors. The acceptance or rejection of jet configurations is realised according to this 
reweighting. Then, each hard parton of the reweighted matrix-element final state undergoes vetoed 
parton showering, i.e. starting from the scale where this parton appeared first, any new emission that 
would give rise to an extra jet is vetoed. In this way, the dependence on the, in principle, arbitrary 
separation scale Qcut regularising the matrix elements is mostly eliminated and the accuracy of the 
parton shower is preserved. However, the cancellations are exact only up to at most next-to-leading 
logarithmic order. This leaves some unavoidable residual dependence, which can be used to tune the 
procedure and obtain optimal agreement with data. 

In Sherpa the CKKW merging of matrix elements and parton showers is accomplished as follows: 

1. All cross sections Ufc for processes with k — 0,1,..., N extra partons are calculated with the 
constraint that the matrix-element final states pass the jet criteria. They are determined by a kx 
measure and the minimal distance is set by the actual merging scale Qcut- Beyond that, Qcut also 
acts as a regulator setting the factorisation (PDF) as well as the renormalisation (as) scales of 
the matrix-element calculations. The kr measure used for jet identification in electron-positron 
collisions can be written as [97| 



and quantifies the fc^ distance between the final-state particles i and j. The jet cuts are satisfied if 
Qij > Qcut- For hadron-hadron collisions, a fc-r scheme is employed, which defines two final-state 
particles to belong to two different jets, if their relative transverse momentum squared, defined as 



is larger than the critical value Q^uf addition, the transverse momentum of each jet has to be 
larger than the merging scale Qcut- The magnitude D, which is of order 1, is a parameter of the 
jet algorithm, see [98j . 




(21) 




(,) (,) 1 2 cosh(y(^) - ) - cos((/)(') - </)(^") ) 



(22) 
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2. Processes of fixed parton multiplicity are chosen with probability Cfc/^^pCr;. The event's hard 
process is picked from the list of partonic processes having the desired multiplicity and according 
to their particular cross-section contributions. All particle momenta are distributed respecting the 
correlations encoded in the matrix elements. Merged samples therefore fully include lepton-jet and 
jet-jet correlations up to N extra jets. 

3. The parton configuration of the matrix element has to be analysed to eventually accomplish the 
reweighting. The partons are clustered backwards according to the same kx jet clustering algo- 
rithm used for the regularisation of the final-state phase space of the matrix elements (cf. step[T]). 
The clustering is guided by the physically allowed parton combinations and automatically yields 
the nodal kx values Qij of each parton emission. It is stopped after a 2 ^ 2 configuration (a 
core process) has been identified. In fact, this backward clustering constructs a limit of leading 
logarithmic accuracy of the full radiation pattern, i.e. determines a possible parton-shower history. 
The sequence of clusterings can thus be taken as a pseudo shower configuration, off which the 
event's evolution will be properly continued by parton showering. 

4. The reweighting proceeds according to the reconstructed shower history. For the strong-coupling 
weight, the identified nodal kx values are taken as scales in the strong-coupling constants and 
replace the predefined choice of the initial generation. The Sudakov weight attached to the matrix 
elements accounts for having no further radiation resolvable at Qcut- The NLL Sudakov form 
factors employed, cf. are defined by 



Ag(Q,Qo) = cxp}- J dqTqiQ,q)[ , 
i Qo ) 

AgiQ,Qo) = cxpl- J dq[Tg{Q,q) + Tfiq)]\ 



(23) 



(24) 



where ^q,g.f are the integrated splitting functions for q — > qg, g gg and g — > qq {ff), given by 



TT q \ q 12 

Tfiq) = (27) 
Stt q 

respectively. T{Q,q) is cut off at zero, such that Aq g{Q,Qo) retains its probability interpretation 
for having no emission resolvable at scale Qq during the evolution from Q to Qo- Hence, A- factors 
are used to reweight in accordance to the appearance of external parton lines. The ratio of two 
Sudakov form factors A{Q , Qq) / A{q, Qq) accounts for the probability of having no emission resolv- 
able at Qo during the evolution from Q to q. Ratio factors thus are employed for the reweighting 
of internal parton lines. In both cases the lower limit is taken to be Qo = Qcut or Qo = D Qcut for 
partons that are clustered to a beam or to another final-state parton, respectively. 

After the reweighting of the matrix- element configuration, the parton shower is invoked starting 
from the pseudo history constructed before. Each parton continues its evolution at the scale where 
it was produced. For the virtuality-ordered parton shower of APACIC++, this scale is given by the 
invariant mass of the mother parton belonging to the QCD splitting, through which the considered 
parton has been initially formed. Recall that this QCD splitting is identified by the "cluster- 
backwards" procedure. For the cases where the considered leg originates from the core process, the 
(hard) scale still needs to be determined. For example, all four partons resulting from a clustering 
that terminated in a pure QCD 2 — s- 2 process will commence their evolution at the corresponding 
hard QCD scale. 



20 




6. In all circumstances parton-shower radiation is subject to the condition that no extra jet is pro- 
duced. Any emission that turns out to be harder than the separation cut Qcut is vetoed. The 
exception to this veto - called highest-multiplicity treatment - is for matrix-element configurations 
with the maximal number TV of extra partons. These cases require the parton shower to cover the 
phase space for more jets than those produced by the matrix elements. To obtain an inclusive N-jet 
prediction, the veto therefore is on parton emissions at scales harder than the softest clustering 
scale, Qsoftest ^ Qcut- Of course, correlations including the N + 1st jet are such only approximately 
taken into account. 

As a final remark it is worthwhile to note that the procedure described here has been fully automated, 
i.e. in Sherpa the merging of the different jet multiplicities is handled on the fly. 

5.2 Applications 
5.2.1 e+e hadrons 

To exemplify the working of the algorithm described above one may consider the simple example of jet 
production in electron-positron annihilation. 

The lowest-order matrix elements account for the k — 2 processes e+e^ qq. They all have the same 
topology and the corresponding Sudakov weights read 

W2 = [A^iE,^,Qcnt)f ■ (28) 

The vetoed parton showers for both the quark and the antiquark start at the scale equal to the centre- 
of-mass energy Ecm of the e+e^ annihilation. A veto is applied for emissions above Qcut- There are two 
topologies for the k — 3 subprocesses, i.e. for e'^e~ — > qqg, corresponding to the emission of a gluon off 
the quark or antiquark, respectively. The related Sudakov weights read 

= A,(gcm,Qcut) ^^^f'^^Q^-'^ A,(g,geut) Ag(g,Qcut) 
^qCl, Vcutj 

= [A,(£;cm,Qcut)]'Ag(g,Qcut), (29) 

where q denotes the scale at which the gluon has been resolved. Again the parton showers for the quark 
and antiquark start at Ecm- The shower for the gluon is initiated at scale q. All showers are subject to 
a veto on Qcut • Note that in the case N — k — 3 the weight reduces to 

W3 = [A,iE,n.,q)f (30) 
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Fig. 9 



Pt [GeV] 



Sherpa CKKW shape predictions for two different merging-scale choices and various pr spectra of 
the vector boson in ^+jets production at Tevatron Run II. The left part shows the distribution 
for the fuU pr range as measured by D0, whereas the other parts focus on the peak region. The 
rightmost plot has been obtained by requiring forward-rapidity lepton pairs. Data are taken from ^5] 
and the shaded bands visualize the sum of statistical and systematic uncertainties on the data. 



and shower vetoes are performed w.r.t. q. 

The quality of the procedure can properly be studied by means of the differential jet rates, which show 
the distributions of y = Q"^ I S values at which an n + 1 jet event is merged into an n jet event according 
to the jet clustering scheme employed by the algorithm. In particular they allow a thorough inspection 
of the transition region y « ?/cut- The merging procedure provides a good description of both the 
hard and the soft part of these distributions. Large values of y indicate the region of hard emissions, 
where the matrix elements improve over any given parton-shower estimate. For decreasing values of y, 
the real-emission matrix elements diverge and descriptions based on them become unreliable. Instead 
kinematically enhanced logarithms need to be resummed to all orders as accomplished by parton showers. 
Shower-level results as given by the procedure for the Durham differential n — > n+ \ jet rates are depicted 
in Fig. [5] for n = 2, 3. They clearly show the features described above: predictions for different values 
of ycut vary at a level of 15% only. The transitions between shower and matrix-element domains for the 
various ycut choices are smooth. The plots emphasise that the merging procedure works well and at the 
level of accuracy it claims to provide for the parton-shower tree-level matrix-element matching. 

5.2.2 W/Z+iets 

Weak boson production in association with jets is one of the most prominent backgrounds in many 
new-physics searches with signals of the type multijets plus missing transverse energy plus leptons. 
The validation and tuning of Monte Carlo tools against newest Tevatron data on weak boson plus jets 
production therefore has to be considered a crucial prerequisite for successfully applying the available 
tools in LHC analyses. The extrapolation to higher energies can only be more meaningful once it 
starts off a solid and well understood test scenario. With more and more Tevatron data coming in, its 
discriminating power gradually increases and data can thus be used to refine the multijet-production 
algorithms, even rule out choices or falsify assumptions that have been made in constructing them. 
Here, new comparisons are presented using very recent experimental W/Z-\-]ets data taken during Run 
II at the Fermilab Tevatron and published by the D0 and CDF collaborations in 99 and [ ^llOOj . 
respectively. All CKKW predictions in this subsection have been obtained from inclusive samples where 
matrix elements including up to three extra (light-flavour) partons were merged with the parton showers. 
The PDF set used was CTEQ6L. Note that effects of the underlying event and hadronisation have not 
been taken into account for this study. 

Figure [H] shows transverse momentum distributions of the Drell-Yan lepton pair. On a rather inclusive 
level these spectra allow for insight to the nature of the occuring QCD radiation, which the boson is 
finally recoiling against. The data are taken from a D0 shape measurement of Z boson events with 
mass 40 < M < 200 GeV and integrated luminosity of almost 1 fb^^ [99J. The overall agreement is 
satisfactory. In the left panel of Fig. [5] the full px range measured up to 260 GeV is depicted. The 
two CKKW predictions differing in their choice of the merging scale Qcut tend to undershoot the data 
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Inclusive n jet cross sections compared to CDF Tevatron Run II data for W^+jets production [92] 
(left panel) and Z+jets production |100] (right panel). The shaded bands represent the full error 
on the data, except for the luminosity error being not included for VK+jets. In both cases the 
blue solid line labels SHERPA CKKW predictions obtained with the CTEQ6L PDF set; for the 
red dashed curve, the CTEQ6M set was used. The predictions for W and Z+jets are normalized 
to the inclusive n — and n = 1 cross sections, resulting in constant if -factors of 1.33 (1.14 
for the dashed curve) and 1.71, respectively. For cuts and details of both analyses, the reader 
is referred to the respective publication. 



between 30 < pr < 80 GeV. This may be improved further by tuning the merging scale and fully 
including underlying-event plus hadronisation effects. The center and right panels respectively show 
the region of soft pt without and with the additional constraint of reconstructing the rapidity of the 
Drell-Yan pair to be larger than 2. For both plots, the results in the very first bin may be corrected by 
tuning the intrinsic fcr smearing parameters. The good agreement, in particular for forward lepton-pairs 
events, confirms that resummation effects are reasonably accounted for. 

The TeV centre-of-mass energies of the Tevatron/LHC provide sufficiently /enormously large phase space 
to be filled by relatively hard QCD radiation giving rise to (hard) jets accompanying the vector boson. 
Such topologies cannot be generated by parton showers. Figure [TU] displays inclusive cross sections 
for various jet multiplicities. These cross sections are very sensitive to the correct description of QCD 
multijet final states and hence provide excellent probes for Sherpa's merging approach. Figure fTTl shows 
the transverse momentum spectra of jets produced in association with the Drell-Yan lepton pairs. These 
observables allow a more detailed study of the QCD structure of the events. In both figures CKKW 
predictions have been compared to CDF data pubfished in [Hlj for VF-|-jets and |100j for Z/7*+jets 
production. The predictions are normalized by a constant X-factor to the lowest available n jet 
cross sections, n = hq. The description of the n > no cross sections then is very satisfactory and the 
predictions are low for the highest jet only. Merging in matrix elements with even more extra partons 
will certainly improve on this situation. The same constant ii'-factors were employed to generate the 
jet Pt spectra. As a result rate and shape of all pr distributions is in quite remarkable agreement with 
the data. There is a tendency of overestimating the hardness of the leading jet, yet the predictions are 
consistent with data. For the case of W^-|-jets production, a second CKKW curve was added obtained 
by employing the CTEQ6M PDF set. This led to a smaller X- factor and softer spectra in general. It 
also provides an estimate of the uncertainty that comes along with the use of different PDF sets. 
A clear deficiency of the parton shower is the limitation in correctly accounting for boson-jet (lepton- 
jet) and jet-jet correlations. This drawback can be overcome by using the CKKW merging approach. 
Examples for Z/'f*+X production at the LHC are given in Fig. [T^l The left panel depicts the pseudo- 
rapidity difference between the lepton pair arising from the decay of the vector boson and the hardest jet 

^The new data has been used to further validate and tune SHERPA by revising and finally changing some of the default 
choices of the implemented merging procedure; the new settings will be made available to the public. 
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Fig. 11 Jet pt spectra compared to CDF Tevatron Run II data for W^+jets production [92] (top row) 
and Z+jets production [lUOj (bottom row) . Labelling, _ft'-factors and cuts are the same as given 
in the caption of Fig. [lOl Note that the top row gives inclusive jet pr spectra of the nth jet; 
the bottom row depicts all-jet pr spectra for events with at least n jets, i.e. a 3 jet event gives 
3 entries to all distributions up to n = 3. 



accompanying the boson. The prediction obtained by solely adding parton showering to the Drell-Yan 
lepton-pair production (here given by the APACIC++ shower) shows a suppression at low |A?7|. The 
correct correlation is encoded in the QCD real-emission matrix elements for Z/j* production and can 
hence be accounted for by the CKKW approach. As a consequence the suppression has disappeared in 
the CKKW merging prediction. The distributions of the spatial separation between the second and third 
hardest jet are given in the right panel of Fig. 1121 The pure shower approach predicts separations that 
are too wide. This again is corrected for in the CKKW procedure by including matrix elements here up 
to three extra parton emissions. 

As demonstrated by these examples CKKW predictions lead to considerable improvements where the 
application of parton showers merely gives approximations or is even misleading. An improved ISR 
description as well as a correct treatment of correlations involving jets can only be achieved by including 
the corresponding matrix elements for extra QCD radiation. Merging algorithms are therefore an ideal 
approach to provide inclusive W/Z+jets samples up to a certain number of jets including effects of hard 
ISR and jet correlations. The comparison to recent Tevatron data proves the relevance of Sherpa's 
merging approach in describing multijet events realistically. 

5.2.3 tt+jets 

It is worthwhile to study the implications of the CKKW formalism for relevant processes at future 
colliders. An example for such a process, which plays a significant role, both as a signal for a bet- 
ter measurement of Standard Model parameters and as a background to new physics searches, is top 
quark pair production. In this publication it serves as an example for the merging of matrix elements 
and parton showers in processes with strongly interacting heavy resonances. Properly implemented real 
next-to-leading order corrections can be crucial for this kind of interaction, as pointed out e.g. in [70j . 
The combination of matrix elements and parton showers in this respect is based on using a (poten- 
tially Breit-Wigner improved) narrow width approximation to compute the corresponding hard matrix 
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Left panel: | Ar;| distributions between the vector boson and the first jet in Z/7*+jets production 
at the LHC. Right panel: Spatial separation between the second and third hardest jet, Ai?23, 
in Z/7*+jets production at the LHC. Dark solid curves represent CKKW predictions including 
matrix elements up to three extra partons. The bright dashed curves are pure shower predictions 
generated by APACIC++. Jets are defined by the Run II fcx-algorithm with parameter D = 0.4 
[98], and the cuts are p^^* > 20 GeV and |r;j'='| < 4.5. 



elements, cf. Sec. |31 Through the identification of the decaying intermediate particle, a parton shower 
can independently be assigned to its production and decay. Upon application of the CKKW merging 
algorithm and when integrating over the phase space of the outgoing particles, jet measures Q between 
strongly interacting particles must be larger than a critical value Qcut- Due to the factorised structure 
of the process, this critical value might be chosen separately per subprocess, i.e. there may be different 
values Qprod = Qcut(prod) and Qdcc,i = Qcut(deCj). 

Technically the merging proceeds as follows: Once a particular momentum configuration is chosen for 
the hard matrix element, the parton-shower history is identified as in the standard merging prescription. 
Since the full amplitude factorises over time-like propagators, only particles belonging to the same sub- 
amplitudes ^prod or Adcc.i Can be combined. When reweighting the hard matrix element with Sudakov 
form factors, the potentially different cutoff values Qcut must be employed. Sudakov reweighting for 
the resonant intermediate particles takes places according to the original prescription. The jet veto is 
applied separately within each part of the parton showers related to a different subamplitude, employ- 
ing the corresponding veto scale Qcut- The highest multiplicity treatment is applied separately in each 
subamplitude. To summarise, the above prescription translates into the CKKW method being applied 
separately and completely independent for each part of the process under consideration. 
For this procedure, reweighting with analytic Sudakov form factors for heavy quarks is necessary. It must 
therefore be assured that those match the respective distributions generated by the parton shower. A 
corresponding comparison is shown in Fig. [131 Predictions from APACIC++ and an analytic calculation 
performed in [lOlj are displayed for the Durham /cT-jet rates at a fictious electron positron collider 
operating at 500 GeV. 

Figure [TJI shows the impact of the description of additional hard radiation through appropriate matrix 
elements in the case of tt-production at the LHC. Subsamples from a given matrix element configuration 
are displayed in colour with the colour code explained above the figure. The following notation is 
employed 

• The first number denotes the additional jet multiplicity in the tt production process. 

• The second number denotes the additional jet multiplicity in the t and t decay process. 

It can be observed that the transverse momentum spectrum of the produced heavy quark pair is signifi- 
cantly enhanced at high values when employing the CKKW prescription. This is easily explained by the 
failure of the parton-shower approach to correctly describe hard initial-state radiation. Also, it is clearly 
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Fig. 13 Comparison of the Durham 2- and 3-jet rates for massive quarks, calculated according to |101] . 

with results from APACIC++. Good agreement is found between the parton-shower result (solid 
histogram) and the analytical calculation (dashed line). For reference, results with massless quarks 
are also shown. 



seen that subsamples including an additional hard jet in the production process (1-0 and 1-1) deform 
the pseudorapidity spectrum of the first additional hard jet (identified through fcx-clustering and Monte 
Carlo truth based ^-tagging). The pure parton-shower approach predicts a dip at central rapidity, while 
the merged sample does not show this feature. 

6 Multiple parton interactions 

The hardest (primary) partonic interaction in hadronic collisions may be accompanied by softer (sec- 
ondary) ones, which may, however, still fall in the realm of perturbative QCD. This is a valid assumption 
and clear evidence has come from experimental studies carried out by the CDF collaboration at the Fer- 
milab Tevatron during Run I and Run II [102pi03j . It was shown that a correct description at the hadron 
level of particle multiplicities and jet activities can only be achieved by Monte Carlo event generators 
that incorporate a model for multiple parton scatterings. 

The basic idea of such a model is to postulate the probability distribution for the occurrence of multiple 
scatterings. According to the model of 3S1 this is given as a function of the non-diffractive cross sec- 
tion fiND, the hard perturbative cross section (Thard and the outgoing transverse momentimr p± in the 
scattering by 

Pip^,b) = /J(6)^^, (31) 

where the prefactors fc (normalisation) and f{b) (proton shape function) incorporate an additional 
impact-parameter dependence of the distribution on an event-by-event basis. In the coUinear factorisation 
approach a phase-space cut on the hard matrix elements, e.g. a minimum jet resolution cut Qcut, bas to 
be introduced to obtain a well-defined differential cross section in perturbation theory. Then dcrhard/dpj_ 
receives its dominant contributions from 2^2 processes and, hence, the definition of p± in Eq. (|3ip is 
unambiguous. Moreover it suffices to consider hard 2 — > 2 QCD processes only. 

The multiple-interactions model of Sherpa is implemented in the module AMISIC++, which generates 
multiple scatterings according to the basic ideas listed above. There are however important details 
where the Sherpa approach deviates from the original formalism presented in |35j . Firstly, all secondary 
interactions are undergoing parton-shower corrections and, therefore, AMISIC++ has been interfaced to the 
parton-shower module APACIC++. A second class of extensions concerns the combination of the CKKW 
merging approach with the modelling of the multiple interactions, which essentially are independent of 
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Fig. 14 Transverse momentum spectrum of the tt-paiT (left panel) and rapidity spectrum of the first extra 
jet (right panel) in pp ^ tt —> W^W~bb + jets events at -^i = 14 TeV. Shown is a comparison 
between the pure parton-shower result (no int. PS), parton shower with radiation off intermediate 
top quarks (Apacic++) and a CKKW-merged sample. Contributions from various matrix element 
configurations are highlighted in colour for the merged sample, with the colour code indicating the 
additional jet multiplicity in the production ® decay process. Up to one extra jet has been simulated 
through matrix elements in the production and each decay. 



the treatment of the hardest scattering in the collision. It is vital that the parton showers related to 
secondary interactions respect the initial p±^ distribution of the hard scatterings. In particular, partons 
with a fcr-separation from other partons larger than p±^ of the respective (secondary) interaction must 
not be radiated. The appropriate way to incorporate this constraint into the parton-shower formalism 
is in fact identical to the realisation of the highest-multiplicity treatment in the CKKW approach. The 
corresponding algorithm then works as follows: 

1. Create a kinematic configuration of the final-state particles according to the hard matrix element 
of the primary interaction and run its parton-shower evolution according to the CKKW formalism 
(cf. SecO. 

2. Employ a /cT-type algorithm in the £^-scheme to cluster the complete final state of the previous 
scattering into a 2 ~> 2 process, similar to how this is done to define starting conditions for showers 
in the CKKW approach. Set the starting scale for multiple interaction evolution to the maximum 
of 

(a) the largest relative transverse momentum between QCD partons that have been combined, or 

(b) the transverse momentum in the 2 — > 2 process, if this process contains a QCD parton in the 
final state. 

3. Select the p± of the next secondary interaction according to Eq. (|3ip . If this happens for the first 
time in the event, select the impact parameter b of the collision as described in [35j . Create the 
kinematic configuration of the secondary interaction. 

4. Set the jet veto scale of the parton shower to the transverse momentum selected in step [3l Start 
the parton shower of the secondary interaction at the QCD scale 

2 2stu 

5. Return to stepO 
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Fig. 15 The left panel displays the average scalar pr sum versus pt of the leading charged particle jet while 
the right panel displays the average charged particle multiplicity versus Acp with respect to the 
leading charged particle jet for pt > 5 GeV of this jet. In both cases Sherpa results are contrasted 
with data from [103) . 



The above algorithm works for pure QCD hard matrix elements as well as for electroweak processes being 
the primary hard scattering in the event. In the QCD case the selected starting scale for the determination 
of the first additional interaction reduces to the transverse momentum in Eq. ([?T|) and is thus equal to 
the original ordering parameter. For electroweak core processes, e.g. W or Z boson production, there 
is no such unique identification. However, it is required that the multiple scatterings in the underlying 
event do not spoil the jet topologies generated by the primary hard scattering. Regarding the electroweak 
bosons as being radiated off QCD partons during the parton shower evolution of a hard QCD event, it 
is appropriate to reinterpret this hard matrix element as a QCD+EW process, whereof the simplest is 
an EW+l-jet process. In this way all primary jet topologies can be preserved. 

An appropriate method to include an impact-parameter dependence in the multiple-interactions model 
is described in [31]. In AMISIC++ the default (anti-)proton shape in position space is chosen to follow a 
Gaussian matter distribution. AmisiC-i~-i- has primarily been validated using the data presented in |103) . 
These tests revealed that effects owing to the variation of the particular shape function are of minor 
importance for the description of the underlying event. It might however turn out that the shape 
function becomes more important when facing more precise data. 

Finally, it has to be stressed that AMISIC++ relies on complete factorisation of the hard (primary) process 
and the additional scatterings in the underlying event. Corrections, which arise from multi-parton PDFs 
as discussed in |104j . are not taken into account. 

In Figure [15] predictions from Sherpa are compared with experimental data from 103 . Both observables 
are very sensitive to the correct description of the underlying event. 

7 Parton-to-hadron fragmentation 

After the parton shower has terminated, a configuration of coloured partons at some low scale of the 
order of a few GeV in transverse momentum emerges. These partons, in order to match experiments, 
have to be translated into hadrons. Since there is no first-principles approach yielding quantitative 
results, hadronisation is achieved by phenomenological models only. Usually, they are based on some 
qualitative ideas on how the parton-to-hadron transition proceeds, like, e.g., local parton-hadron duality 
or infrared safety, defining the model's coarse properties. However, many of the important finer details, 
often related to how flavour is created and distributed in the procedure, are entirely undefined and 
subject to phenomenological parameters only. These are essentially free and must be fixed by extensive 
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comparisons with data, as done for instance in [105j . 



7.1 Beam- remnant handling 

In collisions involving hadronic initial states, there is the additional problem that one or more final-state 
particles will be colour-connected to the hadron remnants. The breakup of the initial hadron into partons 
then needs to be modelled suitably. This is achieved in such a way that only a minimal set of particles 
(quarks and diquarks, the latter as carriers of baryon number) is produced in order to reconstruct 
the constituent flavour configuration of the corresponding hadron. Once their colour connections are 
established, the complete set of final-state partons can be hadronised affecting energy flows and other, 
similar properties of the event at hadron level. 

The distribution of colour in the hadron remnants cannot be inferred from first principles. Hence phe- 
nomenological models have to be employed. In Sherpa the default distribution of colour is guided by 
the idea of minimising the relative transverse momentum of colour dipoles spanned by outgoing partons. 
When including multiple parton interactions in the simulation, it is not always possible to accomplish 
free colour selection in the hard process and minimisation of relative transverse momenta simultaneously. 
In such cases the colour configurations of the matrix elements are kept but the configuration of the beam 
remnants is shuffled at random until a suitable solution is found. The shuffling is implemented also as 
an alternative choice if no multiple interactions are present. So far, no significant differences have been 
found in the predictions from the two models on the level of physical observables. 

In addition to the issues related to colour neutralisation with the beam remnants, all shower initiators 
and beam partons obtain a mild fc^-kick according to a Gaussian distribution parametrising the Fermi 
motion of these particles inside the incoming hadron. It should be stressed that within Sherpa mean and 
width of this distribution are rather small (about 0.2 — 0.8 GeV), such that the physical interpretation 
indeed holds. 

7.2 Cluster-hadronisation model 

For a long time, for the hadronisation of partons, SHERPA has relied on the implementation of the Lund 
string model [106) in Pythia, accessible through a corresponding interface. Since version 1.1, SHERPA 
employs its own module AhaDIC jl07j . which implements a cluster-hadronisation model as described 
in |36j and extended by ideas presented in [37]. The basic assumption underlying this class of models is 
local parton-hadron duality, i.e. the idea that quantum numbers on the hadron level follow very closely 
the fiow of quantum numbers on the parton level. In this framework, the mass spectrum of the emerging 
colour-neutral clusters is dominated by typical hadron masses or masses slightly above. It is therefore 
natural to think of them as some kind of "hadron matter" , carrying the fiavour and momentum quantum 
numbers of hadrons. This motivates to translate the light clusters directly into hadrons or, if they are 
too heavy, to treat them like hitherto unknown heavy hadron resonances, which decay further into lighter 
ones. The idea underlying AhaDIC is to take this interpretation very literal, to compose clusters out of all 
possible flavours, including diquarks, and to have a flavour-dependent transition scale between clusters 
and hadrons. This results in solely translating the very light clusters directly into hadrons, whereas 
slightly heavier clusters experience a competition between being either translated into heavy hadrons or 
being decayed into lighter clusters, see below. In addition, for all decays, a QCD-inspired, dipole-like 
kinematics is chosen. 

In more detail, in AhaDIC, the hadronisation of quarks and gluons proceeds as follows: 
First of all, the four-momenta of all partons are boosted and scaled such that all of them reside on 
constituent-mass shells (m„ = = 0.3 GeV, nis = 0.5 GeV, TOc = 1-8 GeV, mt = 5.2 GeV), whereas 
the gluon remains massless in contrast to the cluster model implemented in HerwiG and HERWIG++. In 
addition to the quarks, further states will emerge, diquarks, which effectively carry baryonic quantum 
numbers. Their masses are typically given by the sum of the constituent masses of the two respective 
quarks making them up, i.e. mqqr — ruq + mqi . In the following, quarks and diquarks, will commonly be 
caUed "flavours". 

At the onset of hadronisation, in cluster models the gluons must decay. This is because in all hadroni- 
sation models, all hadrons are made of flavour constituents only. This forced decay is the reason, why 
an effective gluon mass enters as parameter in HERWIG and HERWIG++, where the gluon decays with 
no spectator involved. In AhaDIC, the gluon's decay is facilitated in a dipole framework, and there will 
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always be a colour-connected spectator to compensate the recoil when the massless gluon decays into 
massive flavours. In these non-perturbative, enforced decays, the transverse momentum measure p± is 
bounded from above by a parameter p'^^^, typically of the order of the parton-shower cut-off scale. 
In addition, in AhaDIC a non-perturbative as coupling has been encoded for all decays (gluons and 
clusters). It is chosen such that it agrees with a measurement from the GDH sum rule. The corresponding 
form given in jl08j has been implemented with all parameters. In addition to this default setting, 
constant- value choices are also available. 

For the production of the flavours, the available phase space leads to kinematic enhancements of lighter- 
flavour pairs in dependence on the mass of the dipole (i.e. pair of gluon and spectator). It is supplemented 
by constant weight factors, Vq, also used in cluster decays to hadrons. When all gluons have been decayed, 
colour-connected flavours are merged into colour-neutral clusters. In contrast to the original version of 
the cluster model, also diquarks are allowed as cluster constituents from the beginning such that in every 
stage of the hadronisation baryonic quantum numbers can be produced. 

In order to guarantee the universality of this hadronisation procedure, it is important that the cluster 
mass distribution is almost independent of the centre-of-mass-energy in otherwise identical processes. 
This has been checked, e.g. for the case of e+e" annihilations into light quarks, cf. Fig.[Tl] The structures 
in the peak region appear owing to masses of the constituent quarks. 

At every step, all clusters are checked whether they fall into one or both of two regimes, namely 

• either the regime where the mass of the cluster is of the order of the mass of the heaviest hadron with 
a flavour wave component identical to the flavour content of the cluster, i.e. where Mc < ruh + 5h, 

• or the regime where the mass of the cluster is lighter than the summed mass of the two heaviest 
hadrons it can decay into by adding a flavour-antiflavour pair to its flavour content, i.e. where 
Mc<mh + mh' + Shh' ■ 

In order to steer this classification further, two offsets Sh and Shh' have been introduced. 
If a cluster is outside both regimes, it will decay further into two clusters. This proceeds by the flavour- 
antiflavour pair firstly emitting a gluon with transverse-momentum choices as in the case of gluon decays 
discussed above, before then the heavier of the two colour dipoles, made of one of the two flavours and 
the gluon, decays further, exactly as above. 

If in contrast a cluster falls into one or both transition regimes, probabilities for the correspondingly 
allowed transitions into one or two clusters are calculated. For the direct C ^ H transitions, this 
probability is given by 

= JM-AM^-mir + imhTh)' ^^^^ 

where nih and Th respectively are the mass and width of the hadron in question, and where ^!qq{h) 
denotes the corresponding flavour wavefunction component of the hadron (e.g. '^uu{'^^) = Vv^)- For 
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decays C H1H2, the probablity reads 



167rA/3 



4TOim2 



Af2 



(34) 



The parameter n can be used to change the relative rate of transitions and decays. In the decay weights, 
the exponent a steers the preference of the decays for heavier hadrons. Vq is the probabihty to create 
an additional g-flavour pair in the field of the decaying cluster. In each channel, the weights are further 
modified by including weights for the various hadron multiplets. In addition, there is a singlet suppression 
factor to reduce transition probabilities to singlet-octet mixed states in the meson sector. 
Once a decay into two hadrons is chosen, the decay kinematics has to be fixed. This is achieved by 
giving the hadrons a transverse momentum w.r.t. the original flavour axis in the cluster rest frame. It is 
distributed according to 

p±_ = \p\ sin 6 where sin 9 = . (35) 
The parameter rj is adjustable and ^ is a uniformly chosen random number between and 1. 



8 Decays of unstable particles 

Unstable particles can be produced in various phases of the event generation, for example in the signal 
process, or after fragmentation, or even from decays themselves. 

For unstable particles generated by the hard process, Sherpa provides the possibility to specify the 
decay chain directly (cf. Sec. [3]). Since the complete process is translated into Feynman diagrams, full 
correlations for all propagators are respected. This approach however is not ideal for particles with 
many decay channels as appearing for example in BSM models. In these cases one often wants to decay 
a particle from the signal process inclusively. Currently this is only possible for the r-leptons, which 
owing to their hadronic decay channels are treated by the hadron-decay module (cf. Sec. 18. ip . A more 
general solution for other fundamental particles is being worked on for a future version of Sherpa. 
On the other hand, unstable hadrons, generated e.g. during fragmentation, and the r-lepton can be 
inclusively decayed in cascade-like processes. Features of this implementation shall be described in the 
following. 

A natural starting point for the description of such cascades is the branching ratios collected in the 
PDG tables [109] . and to merely choose according to them the decay channel for individual hadrons. 
This is not always possible, since in many cases the respective branching ratios do not add up to one. 
There are several options to restore a complete decay table. For many of the light hadrons, especially 
in the cases of the pseudoscalar and vector-meson multiplets and the octet and decuplet baryons, it is 
viable to shghtly rearrange the branching ratios, such that they add up to 100%. Introducing new decay 
channels, "guesstimated" from flavour symmetries and phase-space arguments is an option for some of the 
higher, not well-measured resonances made of light quarks (like, e.g., the tensor mesons). But especially 
for hadrons involving heavy quarks, the sheer number of potentially open decay channels forbids such a 
treatment. Thus, a mixture of explicit hadron decays and transitions of the heavy-hadron state to quarks 
and gluons, invoking again the hadronisation model, turns out to be a much more powerful strategy. 
However, care has to be taken that the hadronic final state emerging after such a partonic decay has 
not already been covered by an explicit hadron-decay channel coexisting with the partonic mode. In 
the hadron-decay framework of Sherpa this is realised by explicitly vetoing all unwanted hadronic final 
states emerging after the hadronisation of a partonic decay. 

Having selected a specific decay mode, its decay kinematics needs to be modelled. In order to go beyond 
an isotropic distribution of decay products in the phase space, explicit matrix elements are employed. 
They are usually inferred from the spin structure and similar symmetries of the initial and final state, 
giving rise to only a few amplitudes. This simple picture can be further refined by invoking form factors 
(FFs) for certain transitions - notable examples are the weak transitions of heavy quarks of the type 
b c or b typically used in semileptonic decays of heavy hadrons (cf. Sec. I8.2.ip . or FFs emerging 
for hadronic currents in r-decays (cf. Sec. 18. ip . 

The complexity of the emerging hadronic final-states and the multitude of hadron decays make it impos- 
sible to implement and calculate matrix elements and phase space for a full final state consisting of stable 
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particles only. It is necessary to resort to the construction of "chains" of subsequent decays. Naively, it 
looks like a reasonable as well as feasible approximation to deal with each decay of such a chain com- 
pletely independently. Nevertheless, looking a bit more carefully reveals that complete independence is 
not always a valid assumption. Clearly, the spin structure of decaying and intermediate particles poten- 
tially leads to non-trivial correlations among them and other intermediate particles possibly emerging 
in the primary decay. The algorithm presented in [110] has been implemented to take such non-trivial 
correlations into account. To illustrate it, consider the decay B Dp with a subsequent p tttt decay. 
Naively, neglecting form factors, the full matrix element for B — s- Dtttt can be described by 



M 



-■'Y.^P^'b+P^) ■ <(Ap)£.(Ap) • K-P^, 



(36) 



This then allows to write the matrix element for B — )■ Dtttt in factorised form as 
M - ^ (B ^ L>p ; Ap) (p ^ TTTT ; Ap) . 



(37) 



Here the only difference to the decay cascade is that the sum over helicities Ap is not computed indepen- 
dently for the two subprocesses. The spin-correlation algorithm takes this into account. 
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Fig. 17 Effect of spin correlations for the angle between the decay planes in h —> t t'^ —> tt VtTt'^i't 



Spin correlations for r-leptons produced by the hard matrix element are also taken into account when 
decaying the r. An example for this can be seen in Fig. [T7] displaying the angle between the decay planes 
of T~ and T"*" produced from a Higgs boson decay. The solid lines represent the analytical predictions 
taken from [lllj . 

Unstable particles such as the p meson have a finite width, which may be quite large. In order to acco- 
modate for this effect, the invariant mass of the decaying particle can a posteriori, i.e. after the particle 
has been produced on its mass shell, be smeared according to a relativistic Breit-Wigner distribution. 
When doing so, however, kinematic bounds must be respected. It is quite obvious that there is an upper 
bound given by the mass of the decaying particle subtracted by the masses of the other decay products, 
which themselves may have to be smeared. Hence, a strategy is needed concerning the sequence of the 
mass modifications. Apart from the upper limit, there is also a lower limit for the invariant mass of the 
decaying particle, given by the minimal mass of its decay products. This inevitably will vary with the 
decay channel - fixing the mass before fixing the decay channel could therefore lead to unwanted biases, 
modifying the branching ratios in an unacceptable way. In Sherpa, these two problems are solved by 
fixing the decay channels of the secondaries before applying the Breit-Wigner smearing on their masses, 
the latter with rising width of the particles. 

This still leaves some room for improvement. Currently, the mass of any unstable particle is distributed 
according to a superposition of Breit-Wigner distributions with different weights given by the respective 
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branching ratios and different lower cut-offs given by the minimal masses of the respective final state. 
This simple picture may be refined by adding in threshold effects emerging when new decay channels are 
becoming kinematically accessible. This can be done by reweigthing with the "true" invariant mass in 
the phase-space factors. For the dicing of the offshell mass of a particle a decaying in a two-particle 
decay a — > be, the emerging correction weight w reads 

^ d^^^Hm'^,ml,ml) ^ m| / A(m;2, to^, Tng) ^^^^ 

with A(a, 5, c) = (a — 6 — c)^ — 46c. For n-particle decays, one can resort to an approximation using the 
mass of the heaviest decay product for mf, and the sum of the remaining decay-product masses for mc. 
In Sherpa these phase-space correction weights for updated masses are used, but no correction for the 
matrix element is applied. 



8.1 r-lepton decays 

The T-lepton is the only lepton being heavy enough to decay into lighter hadrons. Consequently, it 
provides an excellent laboratory for measuring hadronic currents and, additionally, to test the Fermi 
point- like interaction assumption and the CVC (PC AC) hypothesis |112) . 

The mass of the r-lepton is roughly 1.77 GeV, which does not suffice to produce a charmed hadron while 
decaying. Therefore, its decay products are either Icptons or light hadrons consisting merely of up, down, 
and strange quarks. For both, the matrix element of the decay 

T- (P) Vr (k) + products (39) 

is given by 

M = ^u{k,X,^){r-ri5)u{P,Xr) (40) 
where the current depends on the chosen decay channel. 

The branching ratio of the leptonic decay channels is about 35%. This is a pure weak process, which 
can be calculated analytically and serves as a suitable testbed for the left-handed leptonic currents. 
The more involved hadronic channels with more than one hadron have complicated resonance structures 
owing to intermediate particles with a short life time such as the p meson. The corresponding decays 
are described by form factors, which in Sherpa can be parametrised according to the following models: 

• Kiihn-Santaman'a(KS) [113j parametrisation, 

• Resonance Chiral Theory (R^T) [TT4] . 

The KS model is a phenomenological approach describing resonances through their Breit-Wigner form. 
They are introduced "by hand" into the matrix element. RxT on the other hand is an effective field 
theory: the extrapolation of Chiral Perturbation Theory (%PT) to higher energies. Using this approach, 
the decay matrix elements can be calculated analytically including all occuring resonances. 
Such form factors have been implemented for decays into up to three pseudoscalars of the type tt* , , 
K^, and the two modes producing four pions. As an example, Fig. 1181 shows a comparison of four 
different form- factor models with CLEO data [115j . 



8.2 Hadron decays 

With approximately 200 decay tables consisting of 2500 decay channels the majority of the decaying 
particles are hadrons. For some of them, form factors have been very accurately measured, while for 
others not even the branching ratios are well known. In the following the features of hadron-decay 
simulations in SHERPA will be described thereby focussing on heavy mesons. Details about the complete 
package can be found in [38] . 
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Fig. 19 Invariant mass of the lepton-neutrino pair 
for the decay B nVee. 



8.2.1 Heavy mesons 

Mesons containing one heavy quark (mg ^ Aqcd) as a constituent can appear as 

• D+ {cd) and D- {cd), 

• D° (cu) and 0° (cu), 

• Df (cs) and (cs), 

• B+ (bu) and B~ (bu), 

• 5° (bd) and 5° 

• (bs) and BO (bs), 

• B+ (be) and B^ (be). 

In most weak decays, only the heavy quark undergoes the flavour change and enters into the V — A 
current. The second quark only serves as "spectator" being exchangeable among the light quarks. 
For mesons with heavy spectator quarks, additional decay modes may appear, because the spectator 
quark can decay weakly as well. Differences may also appear in the form factors, such that separate ones 
are implemented for decays of the Be meson. 

Excited mesons, like the B* , do not feature very diversified decay channels, but decay to their appropriate 
ground states emitting an additional pion or photon. No form factors have been implemented for these 
two-body decays, but they are treated according to the respective generic matrix element. 
The main decay modes of the pseudoscalar ground states fall into three categories: leptonic decays, 
semileptonic decays and purely hadronic decays. Leptonic decay channels are easily parametrised using 
decay constants and generic matrix elements. 

For semileptonic decays, just like for r decays, it is useful to simply assume a factorised form of the decay 
matrix elements: in this case the heavy-to-heavy or heavy-to-light transitions can be modelled with form 
factors, which are determined by effective fleld-theory approaches such as Heavy Quark Effective Theory 
(HQET) [116] , quark-model predictions [117] or QCD sum rules |118| . Fig. [12] exemplifies the use of 
different form factor models in the decay B —f ttv^g by comparing to predictions by the EvtGen simulation 
package |119j and data from BABAR |120| . For hadronic decays, again a factorisation between the two 
hadronic currents is assumed, using the heavy currents from the semileptonic decays and combining 
them with a hadronic current taken over from r-decays. In fact, for a certain class of heavy-to-light 
transitions, this factorisation ansatz has been proven to yield results correct up to higher orders |121j . 
and it will be used for these and many other decays in SHERPA. 
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8.2.2 Meson oscillations and CP violation 



In addition to spin correlations and finite- widtli effects, more improvements beyond the decay cascade are 
related to the simulation of mixing phenomena in systems of neutral mesons, most notably i?_B-mixing, 
and the modelling of CP-violating effects. 

In the following the relationship between mass eigenstates A'/^/^ and flavour eigenstates M°/M" of a 
neutral meson M is defined as 

I Ml) = p\ M°) + q\ JP) and | Mh) = p\ A/°> - q\ JP) . (41) 

This leads to a non-trivial time evolution M^^^ysit) for flavour eigenstates produced at t = in fragmen- 
tation and decay processes. Four related effects are implemented in SherpA: 

1. Explicit mixing in the event record can be accomplished by setting the parameters Am ~ mn — niL 
and Ar = Th — ^l of the neutral meson. According to the time evolution of the flavour eigenstates, 
the mesons will then decay as the appropriate flavour. 

2. CP violation "in mixing" can be simulated by setting the value of 7^ 1. 

3. Direct CP violation "in decays" can be accounted for by specifying different decay tables for 
particles and their antiparticles. 

4. CP violation "in the interference of mixing and decays" appears in the decay to final states common 
to meson and anti-meson, because of interference terms between the amplitude with a mixed meson 
and the unmixed amplitude. It leads to a time dependent rate asymmetry 

^'^^ ^ nK^ysit) - /) + r(MpV(t) - /) ■ ^ ' 

Defining the quantity Xf — for a decay to the common final state / through the matrix 

elements M, the asymmetry can be expressed as 



S ■ sin(Amt) — C ■ cos(Ar7it) Ar=o 
cosh(Ar|) - sinh(Ar|) 



Ac Pit) = .A^^r^t^ ^ ^ • sin(Amt) - C • cos(Amt) , (43) 

2ImAf 



where S = , , ,^ f, and C ~ , , , denote the mixing-induced and direct contributions to the CP 
asymmetry, respectively. 

The factors S and C have been measured or calculated for many decay channels of the B meson, 
cf. [109j for a listing of measured values. If AF 7^ is chosen, Sherpa calculates A/ from the given 
values of S and C, 



|A/=^, (44) 

Im(A/) = A_ , (45) 
2ReA 



1-C 



f 



= - Vl^^CP^^, (46) 



1 + IA/P 

and simulates according to the non-simplified expression for Acp{t) given in Eq. (I43p . 

A simulation of CP violation owing to interference effects in the decays of B^B^ pairs created from T(4S') 
decays is presented in Fig. [201 Here, one of the B mesons is randomly selected to decay semileptonically 
at time io- This provides a flavour tag at to for the other B, decaying at ti into J /'^Ks- Since both B 
mesons have been produced coherently in the T(45') decay, the reference time for the evolution of both 
flavour eigenstates is set by the semileptonic decay. Thus the rate asymmetry is governed by the time 
difference At = ti — to, which can be positive or negative. 
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Fig. 20 CP violation in the interference between ^ J/*A's and ^ J/^J^s with AT = 0, S = 0.725 
and C = 0. 



9 Multiple emission of photons 

In the previous sections nearly all possible aspects of event generation were adressed. However, higher 
order QED corrections in the form of soft-photon radiation were neglected. Nonetheless, soft-photon 
radiation occurs at any stage of event generation whenever charged particles are involved. They affect, 
e.g., production cross sections, (partial) decay widths of hadrons and r-leptons, and broaden momentum 
and invariant-mass distributions of final-state particles. Therefore these corrections have to be taken 
into account in a realistic simulation. In Sherpa this is accomplished by the PHOTONS++ module, which 
is based on the approach of Yennie, Frautschi and Suura (YFS) [iU] . 

The corresponding algorithm relies on the idea of resumming leading soft logarithms to all orders, rather 
than focusing on leading coUinear terms, which are taken into account in conventional parton showers 
(cf. Sec. H]). Soft logarithms are largely independent of the inner process characteristics and can be 
calculated from the charges of the external particles and their four-momenta only. The YFS formalism 
allows a systematic improvement of the eikonal approximation order-by-order in the QED coupling 
constant - which has rendered it the standard approach for precision-calculations of QED radiation 
effects [122]. 

At present the PHOTONS++ module is only capable to handle single-particle initial states and, hence, 
can only effect these corrections for particle decays. For the hard process, collinearly enhanced photon 
emission off charged leptons and quarks is accounted for during the parton-shower evolution, cf. Sec. [H 



9.1 The formalism 

The YFS formalism shall be reviewed very briefly to illustrate the mechanisms of the resummed soft- 
photon contributions and the order-by-order improvement for the hard emission region. 
Consider a process at zeroth order in the electromagnetic coupling a, described by the matrix element 
A^Q. The sub- and superscripts denote the number of real photons involved and the order of a, respec- 
tively. The four- momenta of the final-state particles are labelled p/, while the incoming momenta are 
denoted by pi. The fully inclusive cross section for having any number of additional soft photons with 
momenta k can then be written as 



E 

71V—0 



(47) 



where ny and nn respectively count the numbers of virtual and real photons involved. 
The starting point of the YFS algorithm is to approximate the dressed matrix elements by the zeroth- 
order expression times eikonal factors that depend on the external momenta only. The correct result 
can then be restored order-by-order in perturbation theory by supplementing the non-leading, process- 
dependent pieces. In the case of adding just one virtual photon this can be formalised as 



Ml 



aBM"o 



(48) 
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with A/q being the infrared-subtracted matrix element inchiding one extra virtual photon. Note that ac- 
cordingly Mq is finite in the limit of vanishing photon momentum, i.e. k 0. All soft divergences owing 
to the virtual-photon insertion are contained in the integrated eikonal B, which is process-independent 
and universal^. When summing over all subsequent insertions of further virtual photons, thereby itera- 
tively applying Eq. ()48|) . the very appealing all-order expression 
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(50) 



can be derived. This can be generalised to any number of real photons by exploiting the unique charac- 
teristics of abelian QED and the fact that virtual photons inserted in closed charged loops do not produce 

any additional infrared singularity. Hence, all Mn^'^'^"'^ are free of soft singularities from virtual-photon 
exchange but may still exhibit such due to real-photon emission. 

However, YFS were able to show that real-photon emission processes can also be factorised - on the level 
of the squared matrix elements rather than the amplitudes. For a single radiated photon, one obtains 
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(51) 



Here, S{k) is an eikonal factor containing the soft divergence related to real-photon emissior0. /3^]^~''"" 
denotes the complete IR-finite (subtracted) squared matrix element for the basic process dressed with 
nji real and ny virtual photons. With Eq. (|51|1 at hand the perturbative series can be reordered such 
that a complete partition of infrared finite and divergent terms is achieved. Moreover, the singular terms 
can be exponentiated by separating the real-emission phase space into a region Q containing the infrared 
divergence such that (1 — SI) is completely IR finite. Introducing the YFS form factoJl 

(54) 



Y{n) = 2a[B + B{n) 



the cancellation of real and virtual infrared divergences can be made explicit. The rearranged form of 
the exact cross section given in Eq. (|Tfl) reads 

oo nR 

^ - E — /d'fpda>U27r)454(^p^_^p^._^fc)e^'(")n^(fc^)Q(^-^) 

rJ^^- ' <J A 1 



'/3o° + 4^ + E®'^^ 



(55) 



The original perturbative series, organised in powers of the electromagnetic coupling constant e and based 
on amplitudes, has been rearranged as a perturbative series in a based on squared matrix elements, whose 
infrared singularities have been extracted. While the squared matrix elements encode all the process- 
specific information, e.g. spin correlations, interferences, and hard-emission properties, the YFS form 
factor describes the resummed universal soft limit. 



®The universal integrated virtual eikonal is a sum over the eikonals of every QED subdipole, i.e. 
S = ^ Bij . 

i<j 

Details can be found in 139] . 

The universal real eikonal is defined as a sum over subdipoles 

S{k) = Y.S^,{k). 
For details see I39| . 

* The integrated real-emission eikonal - in analogy to the universal virtual eikonal B - is defined as 
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— 5(fc)(l-0(fc,n)) , 



2aB(n) = / ■ 
with ©(A:, Q) = 1 if k ^ Q and zero otherwise. 
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In this apparent form of the aU-orders cross section the leading-order process /3g can be factored out 
enabhng the construction of a Monte Carlo algorithm to correct (arbitrary) cross sections or decay rates 
for real and virtual QED radiation. 

9.2 The algorithm 

In Ref. [3S] a concrete version of such an algorithm is presented. Two basic aspects of this ansatz shall 
briefly be discussed here, the definition of the soft region H. and the reconstruction algorithm for the 
leading-order particles' momenta. 

To simplify calculations the singular phase-space region has been chosen to be bounded by an isotropic 
hypersurface in the rest frame of the multipole, i.e. the ensemble of charged particles involved in the 
considered process. It is specified by an energy cut-off lu. Photons outside are generated exclu- 
sively. Their number, energy and angular distribution are determined by J ^jp>5'(fc)0(fc, the soft 
real-emission eikonal. The higher-order hard-emission corrections are then accomodated by introducing 
them as a relative weight. 

This can be achieved by approximating the real-emission matrix elements in the quasi-collinear limit 
using the QED variant of the spin-dependent dipole splitting functions of [123]. However, the terms 
describing the soft limit need to be subtracted as they are already included in the YFS form factor. This 
yields new, soft-subtracted, dipole splitting functions denoted by g^ . For the collinearly approximated 
hard-emission matrix elements one then obtains 

= -^Y.^^^^^^^^^s^'Vg.yfii- (56) 

However, the /3-[ lack any interference terms between the different amplitudes entering. These can only 
be incorporated using the exact real-emission matrix element. This, however, is process specific and, 
hence, is only available for some special cases (a list of available matrix elements is given in Tab. 
Furthermore, to maintain a process-independent formulation as far as possible, extended, composite 
objects, e.g. hadrons, need to be treated as point-like ones. To achieve a higher accuracy in selected 
channels form-factor models can be used. 
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Tab. 3 List of available generic and specific infrared sub- 
tracted squared real-emission (/3i) and virtual- 
correction (/3o) matrix elements {V - vector, F - 
spin- 1 fermion, 5* - scalar). 



A comparison of the impact of the various possibilities to incorporate hard-emission effects can be found 
in Fig. 1211 which depicts the total amount of energy radiated and the angular radiation spectrum in 
leptonic J/i/i-decays, respectively. Both distributions highlight the importance of properly accounting 
for hard radiation and reveal the shortcomings of the approximated matrix elements when compared 
to the exact result. The approximation overestimates hard radiation at large angles, thus enhances the 
total fraction of events with hard photons close to the kinematic limit at half the mass of the decaying 
particle. This effect originates from the missing interference terms and is the more pronounced the 
higher the mass of the emitting particle. The soft limit of course is the same in all approaches, since it 
is incorporated in the YFS form factor Y{fl) to all orders in a. 

The original particles' momenta need to be adapted to satisfy momentum conservation after the higher- 
order QED corrections have been included. This is achieved by distributing the additional photons' 
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E, [GeV] e 

Fig. 21 Left panel: Total photon energy radiated for J/tp{lS) I'^E^ in the J/ip rest frame. Right 
panel: Angular distribution of photon radiation for J/ip{lS) £^£~ in the t^-l~ rest frame. 
The effects of different types of higher-order hard-emission corrections are shown: no corrections 
(dotted), approximated real-emission matrix-element correction (dashed) and exact real-emission 
matrix-element correction (solid). 



momenta uniformly among all, charged and neutral, particles of the process, supplemented by a common 
rescaling of all final-state three-momenta in the multipole's rest frame. This prescription necessitates a 
change of the initial-state momenta as well. However, since these were already fixed when calculating 
the leading-order process, this change is interpreted as a shift of the multipole's rest frame. Additional 
complications arise for initial states involving more than one particle, however, at present PHOTONS++ 
is limited to particle decays. 

On the other hand, photons emitted inside 17 are treated inclusively and although their number is infinite, 
the sum of their momenta is assumed to have a negligible effect on the overall momentum distribution. 

10 Summary 

In this article the original publication on the proof-of concept version of Sherpa [S], dating back to 
2003, has been updated. Since then, Sherpa has matured considerably, especially in terms of physics 
capabilities, reliability and user support. At the same time, however, the original strategy of a separation 
of physics modules and structures for event generation has been maintained. In fact, the underlying idea 
of the physics being confined to individual modules, steered by corresponding handler classes, which in 
turn are employed by appropriate event phase classes during event generation, has proven to be extremely 
fiexible, robust and versatile. Eventually, this construction paradigm has allowed a comparably quick and 
painless incorporation of many new physics aspects and additional user interfaces. This has helped to 
augment SHERPA with many hitherto missing physics aspects, for example hadronisation, hadron decays 
and the radiation of soft photons, thus replacing interfaces to already existing Fortran routines from other 
codefH. These additions render the program a truly self-contained event generator for lepton-lepton and 
hadron-hadron collisions. In addition, many existing physics modules have been upgraded or improved. 
This ranges from the implementation of some models for new physics over a significantly enhanced and 
validated merging prescription for matrix elements and parton showers to the incorporation of a number 
of interfaces to, e.g., the Les Houches accord files for transmitting MSSM parameters or the LHAPDF 
library of PDFs. In addition, the user-interface has significantly improved, now including, among others, 
various event formats, a large variety of phase-space cuts and the greatly enhanced SHERPA-internal 
analysis tool. It is worth stressing that all these improvements certainly enhance both the reliability and 
the physics abilities of SHERPA. 

In addition to the published features of SHERPA described in this article, many other features have 
established, but are not yet part of a publicly available version. These features include two new QCD 
shower modules, one based on Catani-Seymour dipole factorisation [124] and a dipole shower [l25j. 
In the matrix-element sector, an automated Catani-Seymour dipole subtraction algorithm has been 

^ For the sake of backward-compatibility, these interfaces are kept available. 
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implemented within AMEGIC++ [126] . and a new matrix- element generator for very large multiplicities, 
Comix [S7]) has been constructed. 

These new modules, have not been released yet as part of Sherpa, but are planned to be included in a 
version 1.2, which will be made available in 2009. Their incorporation can be seen as a first step of the 
package being upgraded to next-to-leading order accuracy in the matrix- element sector. In this context, 
the initial construction paradigm of Sherpa again proves to yield a very flexible code, which can easily 
support the newly emerging structures. 
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